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We study the dynamical formation of topological defects in nonequilibrium 
phase transitions of classical and quantum field theory. We examine three model 
systems. 1). The symmetry-breaking phase transition of a quantum scalar field 
in a Friedmann- Robert son- Walker universe is analyzed through a first-principles, 
microscopic approach in which the nonperturbative, nonequilibrium dynamics 
of the two-point function is derived from the two-loop, two-particle-irreducible 
closed-time-path effective action. Identifying signatures of correlated domains in 
the infrared portion of the momentum-space power spectrum we find that the av- 
erage domain size scales as a power-law with the expansion rate of the universe. 
The observed power-law scaling, for both overdamped and underdamped evolu- 
tion, is in good agreement with the predictions of the Kibble-Zurek mechanism of 
defect formation and provides evidence of the freeze-out scenario in the context of 
nonequilibrium quantum field theory. 2). The formation and interaction of topo- 
logical textures is analyzed in the nonequilibrium phase transition of a classical 
0(3) scalar field theory in 2 + 1 dimensions. At early times, we observe power- law 
scaling of the length scales of the texture distribution in good agreement with 
the Kibble-Zurek mechanism. However, the scaling of the texture distribution 
changes at late times. We provide quantitative arguments that by the end of the 
transition the length scales of the texture distribution result from a competition 



between the length scale determined at freeze-out and the ordering dynamics of a 
textured system. Therefore, a quantitative understanding of the defect network 
at the end of the phase transition requires an understanding of both critical dy- 
namics and the interactions among topological defects. 3). We discuss a black 
hole phase transition in semiclassical gravity. We review the thermodynamics 
of a black hole system and determine that the phase transition is entropically 
driven. We introduce a quantum atomic model of the equilibrium black hole 
system and show that the phase transition is realized as the abrupt excitation 
of a high energy state. We investigate the nonequilibrium dynamics of the black 
hole phase transition and explore similar examples from the Kosterlitz-Thouless 
transition in condensed matter to the Hagedorn transition in string theory. 
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Chapter 1 



Introduction 



1.1 Phase transitions and topological defects 

Phase transitions are common in the world around us and in the laboratory. 
For example, a ferromagnet exhibits spontaneous magnetization below the Curie 
temperature and Helium-4 exhibits superfluidity below ~ 2K (at ordinary pres- 
sures). In the energetic environment of the early universe phase transitions played 
a decisive role in shaping the cosmos we observe today. These include transitions 
at the Grand Unification (GUT) scale, Tc ~ 10^^ GeV, and at the electroweak 
scale, Tc ~ 10^ GeV. In addition. Quantum Chromodynamics (QCD) is expected 
to display a color deconfinement and/or chiral phase transition at temperatures 
Tc ~ 100 MeV. An inflationary phase transition, a possible explanation for the 
flatness, horizon and monopole problems of the standard Friedmann-Robertson- 
Walker (FRW) cosmology, also may have occurred at temperatures near the GUT 
scale. At still earlier epochs near the Planck energy, Tc ~ 10^^ GeV, the classical 
properties of spacetime described by general relativity may emerge from a phase 
transition in candidate theories of quantum gravity. Understanding the dynamic 
processes by which phase transitions occur is of universal physical importance. 

Phase transitions are often accompanied by broken symmetry. For exam- 
ple, crystalline ice formed when water freezes does not share the rotational and 
translational symmetry of liquid water. If the broken-symmetry phase also has 
a degenerate ground state, then the phase transition can produce topological de- 
fects. Relics of the high-temperature symmetric phase of the system, topological 
defects are topologically stable field configurations that are locally trapped in 
an excited state. Examples of defects include vortices in both type-II supercon- 
ductors and the superfluids, Helium-3 and Helium-4. Vortex topological defects 
also appear as cosmic strings in various models of the early universe. Topological 
defects are classified by the homotopy groups n„(M) of the vacuum manifold, 
M. In three spatial dimensions simple defects are domain walls if IIq{M) ^ 1, 
strings if Ui{M) ^ 1, monopoles if n2(M) =^ 1 and textures if U3{M) ^ 1 [||. 

In some condensed matter systems, topological defects are not merely formed 
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during the phase transition, they are the root cause of it. In two spatial dimen- 
sions a system with 0(2) symmetry undergoes a Kosterhtz-Thouless (KT) tran- 
sition 01 . The KT transition from a quasi-ordered phase to a disordered phase 
occurs when enough thermal energy exists to unbind vortex-antivortex pairs. A 
similar transition occurs in a system with 0(2) symmetry in three spatial dimen- 
sions p, |]. In this case, the disordered phase coincides with the appearance of 
long vortex loops spanning the size of the system. 

Kibble was the first to show that topological defects are a generic feature 
of nonequilibrium phase transitions In the course of a symmetry-breaking 
phase transition, a field decays from an unstable, symmetry-restored state to a 
stable (degenerate) vacuum. The correlation length limits the length scale on 
which the field can choose the same vacuum state. In the laboratory, if the phase 
transition proceeds slowly, the correlation length is bounded only by the size of the 
system; the field will effectively choose a homogeneous vacuum state. In the early 
universe, however, the correlation length is bounded by the size of the particle 
horizon. With correlated regions limited in size by causality. Kibble argued that 
phase transitions in the early universe necessarily leave both a domain structure 
of vacuum states and a network of topological defects. 

Topological defects formed in the early universe can have profound conse- 
quences on the subsequent spacetime evolution. In the standard FRW cosmology, 
GUT scale monopoles, even if formed with an initial density equal to Kibble's 
lower bound of one defect per horizon volume, contribute more than 10^^ times 
the largest possible density of the universe consistent with observations. The 
overproduction of monopoles is a puzzle within the standard FRW cosmology 
and provided part of the original motivation for an early epoch of inflation P|. 
If monopoles are formed in a phase transition before or during an inflationary 
phase, their number density is safely diluted by the exponentially expanding scale 
factor. However, it has also recently been observed that in the very energetic pro- 
cess of reheating, as coherent oscillations of the infiaton decay into other fields, a 
nonthermal spectrum of topological defects may be produced, thus again raising 
the monopole problem 0. The details of the reheating process and the pos- 
sible formation of a defect distribution at the end of an inflationary phase are 
currently an area of active study (see e.g 0). Other topological defects such 
as cosmic strings are viable candidates for the seeds of structure formation and 
produce a characteristic signature in the angular power spectrum of the cosmic 
microwave background (for a recent, short review see 0). Satellites such as MAP 
and Plank designed specifically to measure fluctuations in the cosmic microwave 
background radiation to high accuracy are scheduled to be launched in the near 
future. The data obtained from these missions is expected to clarify the precise 
role of topological defects in the origin of structure formation in the universe. 

Nonequilibrium phase transitions are observed in many laboratory systems 
and the study of topological defects has helped forge links between the study of 
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high energy processes in the early universe and condensed matter systems |Ty] . In 
hquid crystals, the formation of a network of topological string defects has been 
observed in the rapid cooling of the system through its critical point, demon- 
strating not only the viability of the Kibble mechanism in producing topological 
defects, but also, the scaling of the defect network [[TT], |T2|. Other condensed 
matter systems, like superfluid ^He, have complex symmetry-breaking patterns 
and offer many potential analogs to early universe processes . 

The proliferation of field theories containing topological defects and their im- 
portance in both condensed matter and early universe processes provides strong 
motivation to extract predictions from defect models that can be compared with 
experiments and observations. These predictions generally contain three main 
ingredients: the density of topological defects immediately following the forma- 
tion process, such as a phase transition; the subsequent evolution of the defect 
distribution; the calculation of physical observables, such as cosmic microwave 
background radiation anisotropies and polarization. While both numerical and 
analytical work has been done on the evolution of the defect distribution and the 
extraction of physical predictions, less attention has been focused on quantitative 
predictions of the initial defect density. 



1.2 The Kibble and Zurek mechanisms of defect 
formation 

The first estimate of the initial topological defect density formed in a phase 
transition was made by Kibble in a cosmological context ||^. The basic ingredients 
of the Kibble mechanism are causality and the Ginzburg temperature To- The 
Ginzburg temperature is defined as the temperature at which thermal fluctuations 
contain just enough energy for correlated regions of the field to overcome the 
potential energy barrier between inequivalent vacua, 

fcTG~e(TG)'AF(TG), (1.1) 

where AF is the difference in free energy density between the true and false 
vacua and ^ is the equilibrium correlation length. We count a fluctuation over 
the barrier as a single degree of freedom. In the Kibble mechanism, the length 
scale characterizing the initial defect network is set by the equilibrium correlation 
length of the field, evaluated at To- In a recent series of experiments the Kibble 
mechanism was tested in the laboratory [0, |T5[. The results, while confirming 



the production of defects in a symmetry-breaking phase transition, indicate that 
^(Tg) does not set the characteristic length scale of the initial defect distribu- 
tion. These experiments were suggested by Zurek, who criticized Kibble's use of 
equilibrium arguments and the Ginzburg temperature |TB . 
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Combining equilibrium and nonequilibrium ingredients, Zurek offers a "freeze- 
out" proposal to estimate the initial density of defects. In equilibrium, and near 
the critical temperature, the correlation length and the relaxation time of the 
field grow without bound as 

? = ^o|er^ (1.2) 

r = ro|er^ (1.3) 
where e characterizes the proximity to the critical temperature, 

^-^7f^^ (1-4) 

c 

and jj, and i> are critical exponents appropriate for the theory under consideration. 
In the freeze-out scenario, the dynamics of the phase transition occur through 
the time-dependent linear temperature quench, 

6=-, (1.5) 

so that for t < 0, the temperature of the heat bath is above the critical tempera- 
ture and the critical temperature is reached at t — 0. The field starts in thermal 
equilibrium with the heat bath at a temperature above the critical temperature. 
Initially, as the temperature of the bath is lowered adiabatically, the field re- 
mains in local thermal equilibrium. However, the ability of the field to maintain 
equilibrium depends on the relaxation time, which diverges at the critical point. 
The divergence of the equilibrium relaxation time as the heat bath approaches 
the critical temperature is known as critical slowing down. Critical slowing down 
results from the combination of a finite speed of propagation for perturbations of 
the order parameter and a diverging correlation length. As the correlation length 
diverges, small perturbations of the order parameter {e.g., lowering of the temper- 
ature) take longer to propagate over correlated regions, therefore it takes longer 
to reach equilibrium. As the critical temperature is approached from above there 
comes a time t* during the quench when the time remaining before the transition 
equals the equilibrium relaxation time 

\t*\=T{t*). (1.6) 

Beyond this point the correlation length can no longer adjust to follow the chang- 
ing temperature of the bath. At time t* the dynamics of the correlation length 
freezes. The correlation length remains frozen until a time \t*\ after the critical 
temperature is reached. In the freeze-out scenario the correlation length at the 
freeze-out time t* sets the characteristic length scale for the initial defect network. 
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Evaluating the correlation length at e{t*), we find the frozen correlation length, 
and therefore the initial defect density, scale with the quench rate as 



an ~ tq—- (1-7) 

The power-law scaling of the topological defect density with the quench rate has 
been verified using 1-dimensional, 2- dimensional and 3- dimensional simulations of 



phenomenological time- dependent Landau-Ginzburg equations [|T^, |T8|, |T^. Ex- 
perimental efforts to test the Zurek prediction in condensed matter systems are 
wide-ranging but inconclusive. While nonequilibrium phase transitions in ^He 
appear to show qualitative agreement 15], the experiments to date lack re- 



liable error estimates and a mechanism to vary the quench rate. In addition, 
quenches in ^ife and high-temperature superconductors |2l| do not seem to 



produce defect densities at the level expected from the freeze-out mechanism. 

The Kibble prediction for the initial length scale of the defect distribution is 
inconsistent with experiment, yet the alternative freeze-out proposal raises many 
conceptual questions. The freeze-out scenario relies on phenomenological ideas 
derived from experience with classical Landau-Ginzburg systems. It is not clear 
how much of this picture, if any, is applicable to quantum fields that undergo phase 
transitions in the early universe. Quantum fields bring new conceptual issues. In 
particular, consideration of decoherence is necessary to understand when and 
how the quantum system at the late stages of the transition can be described 
by an ensemble of classical defect configurations. Even in the classical context, 
the use of equilibrium critical scaling in a fully dynamical setting is, at best, an 
approximation. Critical dynamics is generally much richer and less universal than 
equilibrium critical behavior p2|. For a system coupled to a rapidly changing 



environment, the correlation length evolves through the dynamical equations of 
motion of the field. In fact, the relaxation time is not always well-defined. In 
addition, a realistic bath consists of interactions between the system and other 
fields in the universe. Physical interactions between the bath and the system 
produce interesting and nonequilibrium behavior such as noise and dissipation 
P3| . The reduction of the complicated bath-system interaction to one parameter, 
the quench time scale, and a prescribed linear time-dependence in the effective 
mass of the time- dependent Landau Ginzburg equation neglects these potentially 
important processes. 



1.3 Foundational issues 

Despite significant effort, the detailed mechanisms responsible for the formation 
of topological defects remain only partially understood. Throughout this disser- 
tation we explore three interwoven and complementary themes that are relevant 
in a broad theoretical context and specifically to the problem of nonequilibrium 



5 



defect formation. These themes are (i) equihbrium and nonequihbrium, (ii) quan- 
tum and classical and (iii) fundamental and phenomenological. Phase transitions 
and topological defects occur in both classical and quantum systems, under equi- 
librium and nonequihbrium conditions and are depicted by phenomenological and 

fundamental theories. The study of the nonequihbrium formation of topological 
defects therefore provides a broad theoretical vantage point from which to further 
elucidate the connections between these complementary and disparate concepts. 

1.3.1 Equilibrium and Nonequilibrium 

The formation of topological defects is a nonequilibrium process. Yet, both the 
Zurek and Kibble mechanisms rely on aspects of equilibrium physics. This is, at 
best, inconsistent. Certainly, neither mechanism applies in a far-from-equilibrium 
phase transition. If, as it seems, some aspects of the phenomenological freeze-out 
scenario are correct then it is important to grasp why and under what physi- 
cal conditions. Ultimately, a complete understanding requires a first-principles 
theory of the dynamics of defect formation. 

1.3.2 Quantum and Classical 

Phase transitions and the formation of topological defects occur in quantum 
environments such as the hot, early universe. However, both the Kibble and 
Zurek mechanisms provide a description of defect formation based entirely on 
classical reasoning. We expect differences in a fully quantum process. In some 
situations, the differences between a quantum and classical system are relatively 
small. For example, quantum theory can modify critical scaling exponents and 
change the coarsening dynamics from that of a classical system. However, the 
conceptual differences between quantum and classical systems can also be vast: 
even the definition of a topological defect in a quantum field theory with arbitrary 
quantum state is unknown. In addition, the formation of topological defects in a 
system of quantum fields involves the element of decoherence; topological defects 
can only be identified as classical (stochastic) field configurations when a positive- 
definite probability distribution emerges from the density matrix of the quantum 
system. 

1.3.3 Fundamental and Phenomenological 

Topological defects appear as classical field configurations in phenomenological 
Landau- Ginzburg theories and in the classical limit of symmetry-broken quantum 
field theories. The existence and significance of topological defects require that 
we enlarge our notions of fundamental degrees of freedom. In particle physics 
and cosmology we customarily describe many systems by linear perturbation of 
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interacting field theories. Topological defects do not survive this reductionist 
approach; they represent an important class of nonperturbative fluctuations. In 
special situations there are exact dualities between field and defect descriptions. 
More often, the degrees of freedom we render important and fundamental, defect 
or field, string or membrane, metric or black hole depends on the nature of the 
questions we ask. 



1.4 Organization 

This dissertation is organized as follows. In Chapter 2 we review the formal- 
ism of nonequilibrium quantum field theory, focusing on the techniques useful 
in describing the dynamics of a second-order phase transition. We derive the 
quantum effective action, giving a diagrammatic loop expansion for the both the 
IPI and 2PI effective action. We then review the closed-time-path (CTP) for- 
malism, describing in detail the differences from the conventional formulation of 
quantum field theory. We present the CTP effective action and give an explicit 
expression at two-loop order. As an example, we apply the CTP and effective 
action formalism to a simple system, the quantum dynamics of an anharmonic 
oscillator. 

In Chapter 3, which describes original work published in Ref. [^, we apply 
the techniques of nonequilibrium quantum field theory to the problem of defect 
formation during a phase transition in the early universe. We study the nonequi- 
librium dynamics leading to the formation of topological defects in a symmetry- 
breaking phase transition of a quantum scalar field with self-interaction in 
a spatially fiat, radiation-dominated FRW Universe. The quantum field is ini- 
tially in a finite-temperature symmetry-restored state and the phase transition 
develops as the universe expands and cools. We present a first-principles, mi- 
croscopic approach in which the nonperturbative, nonequilibrium dynamics of 
the quantum field is derived from the two-loop, two-particle-irreducible closed- 
time-path effective action. We numerically solve the djTiamical equations for 
the two-point function and we identify signatures of correlated domains in the 
infrared portion of the momentum-space power spectrum. During the phase tran- 
sition, the scale factor grows by a factor a final ~ 13 ainuiai for the fastest quench 
and a final ~ 1-3 ainiuai for the slowest quench. We find that the size of corre- 
lated domains formed after the phase transition scales as a power-law with the 
expansion rate of the universe. We calculate the equilibrium critical exponents 
of the correlation length and relaxation time for this model and show that the 
power-law exponent of the domain size, for both overdamped and underdamped 
evolution (determined by the expansion rate r), is in good agreement with the 
freeze-out scenario proposed by Zurek. We also introduce an analytic dynamical 
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model, valid near the critical point, that exhibits the same power-law scaling of 
the domain size with the quench rate. By incorporating the realistic quench of 
the expanding universe our approach illuminates the dynamical mechanisms im- 
portant for topological defect formation and provides a preliminary step towards 
a complete and rigorous picture of defect formation in a second-order phase tran- 
sition of a quantum field. The observed power-law scaling of the size of correlated 
domains, calculated here in a nonequilibrium quantum field theory context, pro- 
vides preliminary evidence for the freeze-out scenario in three spatial dimensions. 

In Chapter 4, which, in part, describes original work published in Ref. |2^, we 
study the formation and interaction of topological textures in a nonequilibrium 
phase transition of an classical 0(3) scalar field theory in 2 + 1 dimensions. We 
present two models of the nonequilibrium phase transition, a dissipative quench 
and a pressure quench. The analysis of the dissipative quench is complicated 
by the presence of large thermal fluctuations and is not yet complete. In the 
pressure quench, the phase transition is triggered through an external, time- 
dependent effective mass, parameterized by quench timescale r. When measured 
near the end of the transition (($^) = 0.9) the average texture separation Lg^p 
and the average texture width scale as Lgep ~ 7-0-39±o-02 g^j-^^j ^ ^o.46±o.o4^ 
significantly larger than the single power-law freeze ~ ^^'"^^ predicted from the 
Kibble-Zurek mechanism. We show that Kibble-Zurek scaling is recovered at 
very early times but that by the end of the transition Lsep{T) and Lyj{T) result 
instead from a competition between the length scale determined at freeze-out 
and the ordering dynamics of a textured system. We offer a simple proposal 
for the dynamics of these length scales: Lgepit) = ^i(t) + Li(t — t freeze) and 
Lw(t) = ^2{t) + L2{t — tfreeze) whcrc ~ and ^2 ~ are determined by 
the freeze-out mechanism. Li{t) = {^i)^t3 and L2{t) = are dynamical length 
scales previously known from phase ordering dynamics, and tfreeze is the freeze- 
out time. We find that Lg^pit) and Ly^{t) fit closely to the length scales observed 
at the end of the transition and yield a = 0.24 ± 0.02 and (3 = 0.22 ± 0.07, in 
good agreement with the Kibble-Zurek mechanism using the (overdamped) mean 
field critical exponents yU = 2z/ = 1. In the context of phase ordering these results 
suggest that the multiple length scales characteristic of the late-time ordering of 
a textured system derive from the critical dynamics of a single nonequilibrium 
correlation length. In the context of defect formation these results imply that 
significant evolution of the defect network can occur before the end of the phase 
transition. Therefore a quantitative understanding of the defect network at the 
end of the phase transition generally requires an understanding of both critical 
dynamics and the interactions among topological defects. 



In Chapter 5, which describes work in progress |26, U^, GHI, we discus black 



hole phase transitions in semiclassical gravity. We review the thermodynamics of 
the black hole phase transition, focusing on the calculation of the semiclassical 
free energy of a black hole in thermal equilibrium. We determine that the phase 
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transition is entropically driven: the large amount of black hole entropy makes the 
transition to a black hole spacetime inevitable. We introduce a quantum atomic 
model of the equilibrium black hole system and show that the phase transition 
is realized as the abrupt excitation of a high energy state. We then discuss the 
noncquilibrium dynamics of the black hole phase transition. We highlight the 
inadequacy of the homogeneous nucleation theory of first-order phase transitions 
and explore similar examples from condensed matter and string theory. 
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Chapter 2 



Nonequilibrium quantum field theory 



2.1 Introduction 

In this chapter we present a first-principles, dynamical treatment of quantum 
field theory. We review the Closed-Time-Path (CTP) formalism of nonequilib- 



rium quantum field theory ||29|, |3T], |3^, |33|, |3^, |3^, |36|, |33, |3^, |3^, ^ with 



particular emphasis on the techniques useful in describing the dynamics of a 
second-order phase transition. There are a number of rapidly evolving high- 
energy systems in which nonequilibrium processes are important. The description 
of these systems requires a fully dynamical, initial value formulation of quantum 
field theory. Examples include the far-from-equilibrium dynamics of the quark- 
gluon plasma now being probed in heavy-ion collisions and the dynamics of 



inflation (see e.g. [|^, ^ and references therein). Furthermore, in the hot early 
universe it is likely that the broken symmetries of the standard model of particle 
physics were partially or totally restored. As the universe expanded and cooled, 
symmetry-breaking phase transitions occurred through the nonequilibrium dy- 



namics of quantum fields (see e.g [§4 



Solving the exact dynamics of a system consisting of interacting quantum 
fields is currently technically impossible. For example, the Heisenberg equations 
of motion for the operator O of a quantum field theory with Hamiltonian if, 

idtO = [0,H], (2.1) 

represent an infinite number of coupled, differential equations at each spacetime 
point. The additional complexity of quantum theory is evident even in the quan- 
tum mechanics of a single point particle. The evolution of a classical particle is 
described in Newtonian mechanics by a single second-order ordinary differential 
equation while a quantum particle is described by the Schroedinger equation, a 
first-order partial differential equation. However, it may be that we don't need 
all of the information contained in the wave function. In fact, in many situations, 
it is satisfactory to follow only the degrees of freedom contained in the expecta- 
tion values of a few well-chosen operators. In this work we focus on obtaining 
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evolution equations for the mean field and the two-point function of a quantum 
field theory. 



Organization 



The organization of this chapter is as follows. In Sec. we derive the quantum 
effective action, giving a diagrammatic loop expansion for the both the one- 
particle-irreducible (IPI) and two-particle-irreducible (2PI) effective action. In 
Sec. p.3| we review the closed-time-path (CTP) formalism, describing in detail 
the differences from the conventional formulation of quantum field theory. We 
then combine this discussion with the results of the previous section to derive 
the closed-time-path effective action. We also give an explicit expression for the 
two-loop truncation of the 2PI-CTP effective action. In Sec. we apply the 
CTP formalism to a simple system, the quantum dynamics of an anharmonic 



oscillator. A summary is given in Sec. |2]5 



2.2 The quantum effective action 

The quantum effective action is a functional of correlation functions whose vari- 
ation yields quantum-corrected equations of motion. Analogous to a free energy 
in statistical mechanics, the quantum effective action is useful in understanding 
the extent to which quantum fluctuations influence the static and dynamic be- 
havior of correlation functions of a quantum system. In particular, suppose we 
are interested in the mean field, 

0(x) = (O+|$^^(x)|O_), (2.2) 

and the two-point function, 

G{x,x') = (0+|T[$h(x)<I>h(x')]|0-). (2.3) 
In the presence of a source J{x) we define the vacuum persistence amplitude 

(0+|0_)_; = Z[J] = [ D0e^[^['^l+^-^(^')*M (2.4) 



and the connected generating functional W[J] via 

Z[J]=e^^W. (2.5) 

The mean field is given by 

6W[J] 



'^("^ 6Jix) 



(2.6) 



j=o 
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while 



^ 5J{x)5J{x') 



(2.7) 

j=o 



If the interaction couphng of the fields is small it is useful to partition the action 
into free and interacting pieces 

S = Sfree + Sint, (2.8) 

SO that 



Since the exponent in the path integral is now Gaussian, the correlation functions 
derived from Z[J] can be computed as a power series in the interaction coupling. 
Such perturbative expansions have proven quite successful in computing scat- 
tering cross sections for weakly-coupled quantum field theories such as quantum 
electrodynamics. However, a number of interesting phenomena in quantum field 
theory are fundamentally nonperturbative in nature. In a phase transition in 
particular, the old vacuum is nonperturbatively different from the new vacuum. 
To analyze such situations we seek a formalism that is both nonperturbative in 
the coupling of the fields yet tractable in terms of quantum interactions. The 
effective action provides one such avenue. 

2.2.1 The IPI effective action 

We define the one-particle irreducible (IPI) effective action r[0] by means of the 
Legendre transform 

V[(t>]^W[J]- J d'yJ{y)<P{y). (2.10) 

The equation of motion satisfied by the mean field is found by variation of the 
IPI effective action. 

Since 

<f>{x) = {^h{x)) = (2.12) 
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then (in the absence of sources) the mean field is determined as a solution to the 
equation 



6r\ 



6(l){x) 



0. 



(2.13) 



At this stage, we have done nothing more than rewrite the generating functional in 
a form appropriate for generating equations of motion satisfied by the mean field. 
We still have the full complexity of the quantum field theory. The advantage of 
this approach is that the effective action T[(j)] can be computed diagrammatically 
to arbitrary order in a loop expansion. The loop expansion is a semiclassical 
expansion where the formal expansion parameter is Planck's constant h. The 
loop expansion provides a calculational technique to sum an infinite class of 
perturbative Feynman diagrams and thus satisfies our nonperturbative criteria. 
The computation of the effective action in the loop expansion proceeds in two 
stages. First W[J] is evaluated to the desired loop order. Then, in order to 
achieve the Legendre transform, we must solve the functional equation Eq. ( p.l2|) 
for J[(p]. This process is illustrated below. 

Expansion of W[J] in powers of h is the loop expansion. We expand the 
classical action about a background field configuration 0o 



and let the background field satisfy 



SS 



1 ^ 6^S 



-J. 



+ 



<f>0 



The connected generating functional is 



W[J] = S[(f)o]+<f)oJ + -ihlndet {tD~^) + 0{h'^) 



where 



6^S 



(2.14) 



(2.15) 



(2.16) 



(2.17) 



00 



The tricky part in the calculation of the effective action is the Legendre transform. 
Without quantum corrections. 



6W[J] 
6J 



(2.18) 



We are thus motivated to define 



+ 



(2.19) 
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where 



5W[J] 
SJ 



(2.20) 



is the mean field we are seeking and 0i (a functional of 0) is at least of order h. 
Inverting Eq. ( p.l5|) to view J as a function of 0o we have 



r[0] = W[J] - 0J[0] = 5[0 + 0i] + (0 + 0i) J[0 + 0i] + -ihlndet {iD~^) 

^ (2.21) 

Since 0i is at least first order in h we can functionally expand about and keep 
terms only to first order. The final result is 



ri-zoop[0] = 5[0] + -z^lndet {iD 



iD 



-1 



(2.22) 
(2.23) 



The generalization of this result to higher loop orders is algebraically involved 
but straightforward. First derived in |^ the effective action to arbitrary loop 
orders is given by 



r[0] = 5[0] + -z/ilndet {iD-^) + ri[ 



(2.24) 



where ri[0] is the sum of all one-particle irreducible vacuum graphs for a theory 
described by the shifted action 



Sshift[4', 0] = 5'[0 + 0] - S[(f)] - 



5^[0] 



(2.25) 



As a concrete example we consider the case of a scalar field theory with potential 
V[0] = im^0^ + |f0'*. The shifted action is given by 



S 



□ + + ^0^ 



+ 



+ 



A 

4!' 



(2.26) 



ri[0] is shown in Fig. |2.1| , in which lines and vertices are given respectively by 
the propagator and interactions of the shifted action, Eq. ( p.26|) . In this case the 
effect of quantum fiuctuations is to change the effective mass of the particles 
and to add a new cubic interaction vertex with a strength proportional to the 
mean field. 
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+ 



+ 




+ 



Figure 2.1: Diagrammatic expansion of Ti[(p] to 3- loop order. 



2.2.2 The 2PI effective action 

In situations where fluctuations about the mean fleld are important (as for exam- 
ple in second-order phase transitions) the IPI effective action and its diagram- 
matic expansion need to be generalized. Inthe IPI effective action the two-point 
correlation function G{x, x') is slaved to the mean field: G{x, x') itself has no in- 
dependent dynamics. Instead, we wish to consider the case where both (j){x) and 
G{x, x') are treated on equal footing. We seek an effective action whose variation 
yields independent dynamical equations for and G{x,x'). 

The generating functional in the presence of sources coupled to both the mean 
field and the two-point function is 

Z[J, = es^t-^'-^l = f D(/,e*Nl+/'^(f)-^(f)+5//<^(2/)-^(2/.^)<^W]. 

(2.27) 

As before, wc make the Lcgendre transformation to define the two-particle irre- 
ducible (2 PI) effective action. 



r[0,G] = W[J,K]- J 0(y)J(y)-l 
- y lG{y,z)K{z,y). 




(p{y)K{y,z)4){z) 



(2.28) 



The equations of motion for the mean field and the two-point function (in the 
absence of sources) are 

smG] 



5(t>{x) 
5G{x,y) 



= 0, 



0. 



(2.29) 
(2.30) 
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+ 




+ 



+ 



Figure 2.2: Diagrammatic expansion of T2[4>, G] to 3-loop order. 



The calculation of the 2PI effective action in a diagrammatic loop expansion 



proceeds along similar lines as the IPI case. Here we simply quote the result |^ 



r[0,G] = + ^Tr In G-^ + '-^TrD-\^)G + TS,G], 



where 



50(a;)50(y) 



(2.31) 



(2.32) 



r2(0, G) is the sum of two-particle irreducible vacuum graphs for a theory with 
propagator G{x, y) and interactions determined from the cubic and higher terms 
in the shifted action Sshift[(l> + 0], Eq. (|2.26|) . Fig. ^]2| shows T2[4>, G] for a scalar 
field theory with 0^ self-interaction. 



2.3 Closed-Time-Path formulation 

In this section we review the closed-time-path (CTP) formalism of nonequilibrium 
quantum field theory. Our discussion follows closely that of Ref . 
the quantum amphtude in the presence of a source J 



. We consider 



Z[J] = {<P',t'\<p,t)j = {<p',t'\Terjf'''"\^,t) 
where \4>,t) is an eigenstate of the Heisenberg field operator at time t 



(2.33) 



(2.34) 

and T denotes time ordering. With vacuum boundary conditions on the fields 
in the far past and the far future, Z[J] is the generator for all vacuum to vacuum 
transition amplitudes. Explicitly, 



(0+|T[$H(a;i)...<l>H(a:„)]|0_) = -{ih)" 



5J{Xi)...6J{Xn) 



(2.35) 
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The cross section for a scattering process can be derived from Z[J] using the LSZ 
reduction formula (see for example |^). As explained in the beginning of the 
chapter there are many situations in quantum field theory where scattering cross 
sections do not provide an adequate physical description. In particular, if we do 
not know the vacuum in the far future, Z[J] is not useful. 

To seek a generator for causal processes we evolve the vacuum state with two 
different sources 

Z[J+, J_] = (0_|0„)j_,j^ = (0_|f [e-n-°°^^-*-]T[e*/-°°c»^+*H]]o_), 

(2.36) 

where T denotes anti-temporal ordering. This is the generating functional for 
vacuum expectation values. For example, 

(0_|T[<|.^(xO...$h(x„)]|0-) = m \' V (2.37) 

All other expectation values can be constructed by taking suitable combinations 
of derivatives and sources. To obtain a path integral representation of the gener- 
ating functional 2'[J_|_, J_] we insert a complete set of states at an arbitrary time 
in the far future. 



(2.38) 



This form immeadiately suggests the path integral representation. 



(2.39) 



The functional integration is taken over fields with boundary conditions appropri- 
ate for vacuum in the far past but equal and otherwise unspecified in the future, 

0±(t^ -oo) ^e^™*, (2.40) 
</)+(t oo) = (f)_{t ^ oo). 

Due to these boundary conditions the path integral does simply factorize and 
propagators beyond the Feynman propagator need to be considered. We illustrate 
this for the case of a free field theory, 

S[(P] = -\j d^x 0(x)(n + m^)(t){x) (2.41) 

To evaluate the generating functional, Eq. ( 2.39| ), we make a change of field 
variables, 



+ (2.42) 
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where 0c satisfies the classical equations of motion, 

(□ + m2)0^(x) = J±(x), (2.43) 
with the correct boundary conditions, Eq. ( |2.4(j| ). The classical solutions are 

0^(a;) = - / d^x'Apix, x')J+{x') + A+(x, x')J^{x'), (2.44) 
d^x'AD{x,x')J+{x') + A-{x,x')J-{x'), (2.45) 



(p_[x) = - 

where Ap and Ao are Greens functions and A"*" and are homogeneous solu- 
tions whose properties we review below. The action is 

+ <P] = j d'^x^[(l)%x)J{x) - 0(x)(n + m2)0(x)]. (2.46) 

Since we have disentangled the boundary conditions, the Gaussian path integral 
can be done exactly, and the free-field CTP generating functional is 

Zfr.ee[J+, J-] = exp ^ j [J+(x)A^(a: - x')J+{x') + J+{x)A+{x - x')J4x') 

-J^{x)A-{x - x')J+{x') - J4x)Ad{x - x')J^{x')] . 

The Closed- Time-Path boundary conditions require the consideration of four two- 
point functions, the Feynman and Dyson propagators and the positive and nega- 
tive frequency Wightman functions. We review the properties of these functions. 
The Feynman propogator is 

p(''-'') = J (27r)>2 _ ^2 + = J (27r)4 p2 _ + ^2 _ 

Ap{x — x') has simple poles in the complex po plane and the po integration can 
easily be done using contour integration. We find 

Af{x-x') = ^^(t-t') /^!^-^e'[P>^-"'')-"(*-*')l (2.48) 

where w'i = jP + m?. When t ^ —00, Ap ^ e™*. The Dyson propogator is 
simply the Feynman propagator with the poles reversed, 

Ad{x-x') = i9{t-t') I ^ ^-i[p{x-x')~n,(t~t')] ^2.49) 

J (^tt) 2,Wp 

J_ if)(f'^f\\ f ^ +i\pix-x')~w(t~t')\ 
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In addition to these propagators there are the positive and negative frequency 
Wight man functions 



^ip(x-x') 



2ni6{p^ - m^)9{po) 



(2.50) 



and 



A (x — x') 



(27r)3 2wp 
d^p 1 



,ji[w{t—t')—p-(x—x') 



^-i[w(t-t')-p-(x-x')] 



(2.51) 



Since the free field generating functional is exact these four functions are the four 
quantum expectation values, 



iAp{x — X 
-iAfj^x — x' 
-iA'^{x — x' 

iA~~ix — X 



hG++ix,x') = (0_|r[<l>H(x)$^^(x')] |0_), 

hG^4x,x') = (o_|r[<i>j^(x)$//(x')] |o_), 

hG+4x,x') = (0_|$^,(x)$h(x')|0_), 



(2.52) 



(0_|$^^(x')<f//(x)|0_). 



As with the in-out formalism the CTP generating functional of an interacting 
field theory is built from Zfree[J+, J~]- For a scalar field theory with quadratic 
self-interaction, 



Zint[J+, J-] = exp 



-iX 

IT 



d X 



5JAxY 5JJxY 



Zfree[J+, J-]- 



(2.53) 



The diagrammatic rules of the CTP formalism follow from Eqns. ( |2.53| ) and 



2.52| ). There are two external legs, one corresponding to J+ the other to J_ 



and two interaction vertices with coupling ±A. Connecting the external legs 
and the vertices are the four Greens functions, Eq. (|2.52| ). Each diagram in 
the conventional (in-out) formalism of quantum field theory is expanded in CTP 
language by writing all possible (±) combinations of the external legs, vertices 
and Greens functions. 



2.3.1 The 2PI-CTP effective action 

The CTP effective action combines the nonperturbative nature of the quantum 
effective action considered in Sec. p.2| with the dynamical content of the CTP 
formalism. The 2PI-CTP effective action generates real and causal equations of 
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motion for the mean field {^H{x,t)) and the four two-point functions Gab{x,x' 
of the quantum field theory. The 2PI-CTP equations of motion are 

dTcTp[(l), Gab, 



S(j)+{x) 

STcTp[(j),Gab] 



(</,+=0_) 



5Gabix,y) 



0, 

: 0. 



(2.54) 
(2.55) 



The 2PI-CTP effective action is constructed from Eq. (|2.31| ) and Fig. |2]^ by 
expanding each diagram in all possible (±) combinations of vertices and Greens 
functions. For definiteness and since we will use this expression in the following 
chapter we give the analytic expression for a two-loop truncation of the 2PI-CTP 
effective action for a scalar field theory with quartic self-interaction, 



^CTp[<Pa, Gab] 



+ 



S[(t)+]-S[<P^]-jlndetGab 

d X I d X j\ (x, X )(j"a^(x , x) 



where 



iA'''{x,x') 



^2[4^a, Gab], 

S'{S[<P^]-S[<P^]) 



6(f)aix)5(l)bix') 
and a, 6 = {+, — }. The two-loop contribution is 

^2[(pa,Gab] = — d'^x[G1_^_{x,x)-G'i_{x,x 



(2.56) 
(2.57) 



+ 
+ 



4 J 


I d^x 


4 j 


f d^x 


2 j 


f d'^x 



(2.58) 



2.4 An example: The quantum anharmonic os- 
cillator 

The preceding discussion was somewhat formal. To illustrate OTP and effec- 
tive action techniques we derive the one-loop equations of motion for the mean 
position of a quantum anharmonic oscillator with classical action 



(2.59) 
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The IPI-CTP effective action truncated at one loop is 

ih 



T[x+,x_] = S[x+] - S[x^] - —TrlnA"', (2.60) 



where 

6^S 



iA'''' = T^^. (2.61) 



The nonzero components of are 



a++(t, H) = -(^^ + w^ + 3Xxl^ 6{t - t'), (2.62) 



and 



iA--{t,t') = (^^^ + + 3\xl^ 5{t-t'). 



(2.63) 



The equation of motion for the mean position of the oscillator is given by 

'5r[x+, X- 



5xa 



0. (2.64) 



Using 51ny4 = A ^5A^ variation of the log term gives {A ^)^^{t,t)6iXx{t) and 
the equation of motion for the mean position is 

(^^ + ^^ + ^a;' + 3Xh{A-Y+{t, t)^ x{t) = 0. (2.65) 

The inverse of the operator (A~^ )"*"+, 

G++(A-i)++ = 1, (2.66) 

is the Feynman propagator for fluctuations about the mean field x(t), 

hG++{t,t') = (0|T [6xH{t),6xH{t')]\0), (2.67) 

where Sxnit) = XH{t)—x{t) is the fluctuation operator. We can rewrite Eq. ( p.65|) 
as a pair of coupled equations 

+ -"^^ + + 3A;iG'++(t, t)^ x{t) = (2.68) 
+ + '^^') '^^^(^' = -'^(^ - (2-69) 
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As they stand these equations are not in a convenient form for analysis (in par- 
ticular for numerical solution). We expand the Heisenberg fluctuation opera- 
tor Sxnit) in terms of time-dependent complex mode functions f{t) and time- 
independent creation and annihilation operators, 

6xH{t) = f{t)a + r{t)al (2.70) 

The creation and annihilation operators are defined with respect to the initial 
vacuum 

a|0) = 0, (2.71) 
and obey the usual canonical commutation relations 

[a,d^] = l. (2.72) 
Since [x, p] = ih, the mode functions must satisfy the Wronskian condition 

f{t)r{t)-r{t)f{t) = zh. (2.73) 

Using Eq. (|2.70| ) the Feynman propagator is simply expressed in terms of the 
mode functions 

hG++{t,t') = 0{t - t')f{t)r{t') + e{t' - t) fit) fit'). {2.7 A) 

The IPI-CTP equations of motion for the mean field x{t) = {xh) and the fluc- 
tuation mode functions are 

^ + w' + Xx^ + 3X\f{t)\'^x{t)=0 (2.75) 

^,+w' + 3Xx'^ fit) = 

The equations of motion preserve the Wronskian condition, Eq. (|2.73|) so that 
we only require it to be satisfied initially. The nonlinear (ordinary) differential 
equations, Eq. (|2.75|) , are now easy to treat numerically. 



2.5 Summary 

In this chapter, we reviewed the Closed-Time-Path (CTP) formulation of nonequi- 
librium quantum field theory. We introduced the quantum effective action and 
presented a diagrammatic loop expansion of both the one-particle-irreducible 
(IPI) and two-particle-irreducible (2PI) effective action. We then combined the 
CTP formalism with effective action techniques and gave an explicit expression 
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for the 2-loop truncation of the 2PI-CTP effective action, fn Chapter 3 we use 
this expression to derive a dynamical equation for the two-point function of a 
quantum scalar field undergoing a symmetry-breaking phase transition in the 
early universe. Finally, as a simple example, we used the IPI-CTP effective ac- 
tion to derive evolution equations for the mean position of a quantum anharmonic 
oscillator. 
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Chapter 3 



Quantum critical dynamics: Formation of 
domains in the early universe 



3.1 Introduction 

A complete understanding of the physical issues involved in the formation of 
topological defects in a quantum field theoretic phase transition requires a first- 
principle approach to the nonequilibrium dynamics of quantum fields, a realistic 
treatment of the interaction of the quantum field with other fields that consti- 
tute an environment, and the identification of classical defect configurations from 
the quantum field system. In this chapter, which represents work published in 
Ref . , we analyze the nonequilibrium quantum dynamics of the phase transi- 



tion, focusing on the formation of correlated domains of true vacuum. We choose 
as our model a quantum scalar field evolving through a second-order phase tran- 
sition, initiated by the cooling of a radiation-dominated Friedmann-Robertson- 
Walker (FRW) Universe. We leave to future research a thorough analysis of the 
field-bath interaction and the quantum-to-classical transition (for recent, related 
work on decoherence see [iS|). 



While symmetry restoration in finite-temperature quantum field theory has 
been known since the early work of Kirzhnitz and Linde most previous efforts 



have focused on the equilibrium aspects of the transition. Techniques such as 
the finite-temperature effective potential |5^ and the renormalization group 



51, 52] have been developed to deduce equilibrium critical properties such as the 
order of the phase transition and its critical temperature. However, equilibrium 
techniques are inadequate to study the dynamics of the phase transition. The 
use of equations of motion for the mean field derived from the finite-temperature 
effective potential was criticized in [Q. In general the use of equations generated 
from the finite-temperature effective potential in a dynamical setting results in 
unphysical (acausal and complex) solutions. Although equilibrium techniques 
are clearly inappropriate, solving the full equations of motion for an interacting 
quantum field theory is generally impossible, even numerically. As discussed in 
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the previous chapter, this difficuhy is partially overcome by the development 
of approximation schemes which allow for the evolution of a restricted set of 
correlation functions of the quantum field theory. These methods have been 
applied to a number of dynamical problems in quantum field theory including 
the dynamics of second-order phase transitions ^ . 

The problem of defect formation in a nonequilibrium second-order phase tran- 



sition of a quantum field has also recently received attention p^, The 

results are promising but the studies are incomplete. In these previous approaches 
the phase transition is incorporated through an ad hoc time dependence of the ef- 
fective mass of a free field theory: an instability in the theory is induced when the 
mass becomes tachyonic. The use of a prescribed time dependence of the effective 
mass, while providing a convenient analytic model, lacks physical justification. 
The neglect of interactions confines the applicability of these approaches to very 
early times before the field amplitude grows substantially. They are therefore 
unable to account for the back reaction which is necessary to stabilize domain 
growth and shut off the spinodal instabilities of the phase transition. In addition, 
defects formed during the linear stages of the phase transition are transient and 
not likely to survive to late times. 

It is useful to contrast our first-principles approach to the dynamics of a 
phase transition in quantum field theory with an approach common to condensed 
matter. In condensed matter systems it is common to model critical dynamics 
with a classical, phenomenological, time-dependent Landau-Ginzburg equation 
for an order parameter ^P, 

a,v^(x,t) = -r^ + e, (3.1) 

where F[$] is a phenomenological free energy density for the order parameter, F 
is a phenomenological dissipative coefficient and ^ is a stochastic term incorpo- 
rating thermal fluctuations of the environment [^. Even if the order parameter 
is of quantum origin, as e.g. in the phase of the wavefunction for liquid Helium- 
4, the Landau-Ginzburg equation is rarely derived logically from the underlying 
quantum dynamics of the system. In condensed matter systems, to compensate 
for insufficient microscopic information, the order parameter and its equation of 
motion are chosen with great care and physical intuition. In experimentally inac- 
cessible environments, such as the early universe, it is not a priori obvious what 
the order parameter, or its dynamics, should be. In situations where phenomeno- 
logical approaches are inadequate, it is necessary to work with the fundamental 
quantum dynamics of the fields. A first-principles approach to the quantum dy- 
namics of phase transitions avoids ad hoc assumptions about the dynamics of the 
correlation length and the effect of the quench. The system simply evolves under 
the true microscopic equations of motion. We can therefore explore many details 
of critical dynamics that are inaccessible in phenomenological theories. 
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This chapter is organized as follows. In Sec. p.2| we present the model, a 
derivation of the equations of motion for the two-point function of the theory, 
and discuss renormalization, initial conditions and numerical parameters used 
in the numerical simulation. We then provide a dynamical description of the 
phase transition using results of the numerical simulations. Section begins the 
discussion of domains and presents the argument that domains are determined by 
a peak in the Fourier space structure function k'^G{k, t). Section |3]^ discusses the 
power-law scaling of the size of domains with the quench rate. The equilibrium 
critical exponents of the correlation length and relaxation time are calculated in 
both the underdamped and overdamped cases and the power-law scaling of the 
defect density with the quench rate predicted by the freeze-out proposal is shown 
to be in good agreement with the numerical simulations. An analytical model 
valid for slow quenches and near the onset of the instability is introduced and 
the power-law exponent in the analytic model is found to be the same as the 
numerical simulations. Section |3]^ provides a summary and discussion of these 
results and presents possible directions for further research. 



3.2 Nonequilibrium symmetry-breaking through 
correlation dynamics 

We consider a quantum scalar field in a FRW spacetime. The field has the 
symmetry-breaking classical action 

S = J {^^,^8^^ + m2$2 - ^ ^3 2) 

where g is the determinant of the metric of the classical background spacetime. 
We assume that the stress-energy tensor is dominated by other radiation fields 
present in the early Universe. These fields maintain the overall homogeneity and 
isotropy of the universe. Small deviations in homogeneity and isotropy that may 
eventually be responsible for the fluctuations observed in the cosmic microwave 
background radiation are produced in our model by the topological defects of the 
system and appear only at the end of the phase transition. We therefore work in 
the semiclassical test-field approximation (ignoring back reaction of the $ field 
on the spacetime), and we assume that the scale factor has the time-dependence 
of a homogeneous and isotropic, spatially flat, radiation-dominated universe. 



a{t) 



t + T 



(3.3) 



The expansion of the universe and the resulting redshifting of the modes act here 
as a physical quench allowing the dynamics of the phase transition to unfold 
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naturally. This is in distinction to work which uses an instantaneous change in 
the sign of the square of the mass |^ p9 |. 



The nonequilibrium dynamics of the quantum field is derived through the 
the use of the 2PI-CTP effective action, discussed in the previous chapter. The 
equations of motion for the mean field and the two-point function derived from the 
2PI-CTP effective action respect the $ — >■ — $ symmetry of the classical action. 
Since the field starts in a symmetry-restored state above the critical point where 

{^)^nU^al = 0, (3.4) 

the mean field remains identically zero throughout the phase transition; the dy- 
namics of the phase transition unfold through the dynamics of the two-point 
correlation functions. The equations of motion for the two-point functions are 
obtained through the variation of the two-loop truncation of the 2PI-CTP ef- 
fective action presented in Sec. |2.3.1| . The two-loop truncation of the 2PI-CTP 
effective action is equivalent to the time-dependent Hartree-Fock approximation. 
The equation of motion for the Feynman propagator G{x,x') is 

tt:? T + 3-— -m + —G{x, X) G(x, X ) = -ib[x - x). 

\ot'^ aat I ) ^ 



Instead of directly solving Eq. ( p.5|) we derive the homogeneous equation for the 
(complex) mode functions of the quantum field. We define the spatial Fourier 
transform, 

G(x,x') = 1 ^e-*(^^-^^')Gfc(t,t'), (3.6) 

— * 

where k = \k\, and expand the Heisenberg Field operator $j|/(x, t) as 

^h{x, t)= I K-'hm + e-^'-'frna^^ ■ (3.7) 



The creation operators and annihilation operators and at are defined with 
respect to the initial vacuum, 

%H0) = 0, (3.8) 
and obey the commutation relations 

The field-momentum canonical commutation relations and Eq. ( p.9|) imply the 
Wronskian condition on the mode functions, 

mfm-fmm = ^h. (3.10) 
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The momentum space two-point function Gk{t, t') expressed in terms of the mode 
functions is 



hG,{t,t') = te{t-t'){mfi{t'){N, + i) + fi{t)U{t')N,) 

+ te{t' - 1) {fi{t)Ut'){N, + 1) + Ut)fmN,) , (3.11) 

where 

Nk = (ak^ak) (3.12) 



is the number of particles of momentum k in the initial state. Equations ( p.5|) 
and ( |3.11| ) lead to the mode function equation, 

/ d'^ d(t) d k"^ 9 ^A^, ,\ , , , , , 

where 

/d^k 
j^,h{t)m)a,{T) (3.14) 

is the equal-time limit of the two-point correlation function and 

2iV, + 1 ^ a,{T) = coth (^^^) (3.15) 

is a constant factor incorporating the thermal initial conditions with temperature 
T. The effective mass of the system is 

mlf^{t) = -m' + ^G{t,t) (3.16) 
and the initial frequency is 

wl{0) = e + mljf{0), (3.17) 

where mljj-{0) is the initial finite-temperature effective mass in the Hartree-Fock 
approximation. Since we work in a radiation-dominated FRW universe, the scalar 
curvature, R, is zero and the conformal coupling constant niay be ignored. The 
dynamics of the phase transition is analyzed by solving the coupled, nonlinear 
mode function equations, Eq. (|3.13|) . Before proceeding to this discussion we first 
treat the issues of renormalization, initial conditions and the parameters of the 
model. 
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3.2.1 Renormalization 



The equal-time limit of the two-point function, Eq. ( p.l4|) , appearing in the equa- 
tion for the effective mass, Eq. ( p. 161) , is divergent and must be regularized. A 
simple regularization method, amenable to a numerical simulation, is to imple- 
ment an ultraviolet cutoff in physical spatial momentum. After suitable renor- 
malization, the expression for the two-point function must be independent of this 
cutoff. The cutoff dependence of the bare variance is isolated by considering a 
WKB-type solution to the mode function equation. Identifying the second order 
adiabatic mode functions from the WKB solution and in the limit of large cutoff 
A we find 

The removal of these divergent pieces in the effective mass is implemented through 
mass and coupling constant renormalization. The quadratic divergence is asso- 
ciated with the divergence of the propagator self-energy which we cancel with 
renormalization of the bare mass. We remove the remaining logarithmic diver- 
gence by renormalization of the bare coupling constant. Explicitly, 



2 fi'Xl) A^ 2 

lOTT^ a^(t) 

At- 



(3.19) 



= ' , ^ . (3.20) 

l-^ln(A/«:) ^ ' 

The shift in bare parameters is time independent as long as the cutoff A and the 
renormalization scale k are implemented in terms of the physical momentum, 

A = Aoa(t), (3.21) 
K = Koa{t), (3.22) 

where A and k are comoving and Aq and Kq are physical quantities. The renor- 
malized mode function equation is 

where Gs{t,t) is the subtracted two-point function. 

Our renormalization procedure simply removes the (divergent) contribution of 
short-wavelength modes while leaving the long-wavelength modes important to 
the dynamics of phase transition unaffected. In following sections we will drop 
the renormalization subscripts for clarity and it is to be understood that we are 
working with renormalized quantities. 
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3.2.2 Initial conditions 



The quantum field is assumed to be in an initial state of thermal equilibrium. 
In this model with a tachyonic tree-level mass, symmetry is restored by finite 
temperature corrections. The initial effective mass rn^jj- is the solution of the 
(renormalized) equation 



= -mf ^ -j k^dk \ + e{k - k) " 



In the high temperature and small A limit this yields 

AT2 



(3.25) 



<ff = -m' + (3.26) 

a result familiar from finite-temperature field theory 

In an expanding FRW Universe, exact thermal equilibrium will persist only 
for conformally invariant fields. If the expansion rate is small relative to internal 
collisional processes of the field then there is an approximate notion of equilibrium 
5B|, B3. This is evidenced by transforming to conformal time, ?7, defined by 



dt = a{ri)dri, (3.27) 
and performing a mode redefinition 

fkiv) = fkivHv) (3.28) 
The rescaled conformal mode function equation is 

/ ^2 fc\ \ ^ 

[d^2+k' + «'(^)(^' + Y^(^' V))) MV) = 0. (3.29) 

If the Universe is slowly expanding, a WKB-type solution is appropriate and 
a low-adiabaticity truncation of the instantaneous WKB frequency is sufficient. 
The zeroth-adiabatic order solution to equation ( |3.29 ) is given by 

f{r^) = ^=e-'^'^'='^\ (3.30) 
k"^ + a{rifmlff. (3.31) 



w 



This leads to the following initial conditions for the mode functions in cosmic 
time t 

/fc(0) = ^ (3.32) 

a(0)v/2^' ^ ' 

A(0) = [-"^^-'Mmfm- (3.33) 
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Slow expansion assumes that the natural frequency of the kth mode is faster than 
the expansion rate of the Universe, 

WkiO) > ^. (3.34) 
Zt 

The adiabatic equilibrium approximation fails for the lowest k modes. However 
with high temperature initial conditions the low k modes are not a dominant part 
of the spectrum. 



3.2.3 Numerical parameters 



The model described by Eqns. (|3.3| ), ( p.23|) and (|3.24|) is characterized by 6 pa- 



rameters, (m, r. A, T, A, k). In the following we set h = 1. Spacetime scales 
are measured in units of 1/m, effectively setting = 1. The initial rate of 
expansion or the initial Hubble constant is controlled by r, 

H(0) = ^ = ^. (3.35) 

In our simulations the range of the r parameter is 

0.01 <r< 100. (3.36) 

With this range, the scale factor grows during the phase transition by a factor 
afinai ~ 13 ainitial foT the fastest quench and ajinai ~ 1-3 ainuiai for the slowest 
quench. The self-coupling is (relatively) strong, 

A = 0.1. (3.37) 

The initial temperature T is chosen so that the value of the initial effective mass 
mgjj(O) is of order unity. Specifically we choose 

T = 20.0, (3.38) 

so that 

mljfiO) = -1.0 + = 0.607. (3.39) 

The values of the initial effective mass mgjj(O), the coupling A and expansion 
rate r, were chosen so that the simulations of the phase transition completed on 
numerically accessible timescales. Due to the renormalization scheme, both the 
comoving cutoff A and the comoving renormalization scale n increase with the 
scale factor. The initial value 

Ao = 340 (3.40) 
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was chosen as the lowest value that maintained cutoff independence of the ini- 
tial effective mass. Insensitivity to the cut-off in the dynamical simulations was 
verified by doubling Aq and observing no change in the output plots of the time- 
dependent effective mass. The initial value of the renormalization scale 

Ko = 1.0 (3.41) 

was chosen so that the renormalization scale was always above the maximum 
momentum of the unstable modes. The coupled, nonlinear system of mode func- 
tion equations with the given initial conditions was solved numerically using an 
adaptive stepsize, fifth-order Runge-Kutta code. Mode integrals were performed 
using a simple Simpson rule with a uniform momentum binning 

k = nkbin, (3.42) 

hin = (3.43) 

where Lq = 100.0 is the effective size of the system and n is the total number 
of modes. Insensitivity to the momentum binning was verified by reducing kun 
and observing no change in output. The number of modes n varied from 10^ to 
10^. Run times varied from hours to days on a DEC 500 MHz workstation which 
corresponds to dynamical time scales of t = 2 to t = 100. 



3.3 Domain formation 



The results of a typical simulation are shown in Figs. p.l| and p.2| . In Fig. 
we plot the renormalized effective mass as a function of cosmological time for 
quench parameter r = 1. The phase transition begins when ^e// ^^^^ becomes 
negative. When mljjr is negative, modes with physical momentum ^ < rUeff 
have imaginary frequencies and begin to grow. This indicates the onset of the 
spinodal instability which is characteristic of a second-order phase transition. 
In the early stages of the phase transition, the evolution of the effective mass is 
dominated by the redshifting of stable modes and the effective mass decreases. As 
the phase transition proceeds, more modes redshift into the unstable momentum 
band and the amplitude of unstable modes continues to grow. Eventually, the 
redshift of stable modes balances the growth of unstable modes and the effective 
mass increases. As the effective mass passes through zero from below, the field 
reaches the spinodal point. Beyond the spinodal point all modes are stable. As 
the effective mass continues to grow, non-linear thermal and other coUisional 
processes are expected to be important and the Hartree-Fock approximation of 
the dynamics of the two-point correlation function breaks down (see for example 
5^). In Fig.^]2| we show the Fourier space structure factor 

S{k,t) = eGk{t,t) (3.44) 
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Figure 3.1: Plot of the square of the (renormahzed) effective mass M = m\f^ vs. 
cosmological time t for quench parameter t — 1. 



at various times during the phase transition. The bottom curve is a snapshot of 
S{k, t) at time t = 3.9. The middle curve is a snapshot at t = 5.4. The top curve 
is a snapshot at t = 6.9 As argued in the next section, defects and domains can 
be identified in the low-k structure of S{k,t). As the phase transition begins, 
S{k, t) develops a peak at low k. As the phase transition proceeds, this peak 
grows in amplitude and redshifts until at late times it completely dominates the 
infrared portion of the spectrum. 

3.3.1 Defect density 

We have argued that the dynamics of a phase transition of a quantum field may 
be approximated by the dynamics of the two-point function at least for times 
less than the spinodal time. Using the two-loop 2PI-CTP equations of motion 
we obtained a numerical solution for the evolution of the two-point function. To 
observe the formation of correlated domains and topological defects it is necessary 
to identify these structures from the form of the two-point function. 

The identification of topological defects from the underlying quantum dynam- 
ics of the two-point function is a complicated problem. Intuitively, the existence 
of well-defined correlated domains of true vacuum at the completion of the phase 
transition is similar to the domains that form in a condensed matter system like a 
ferromagnet. However, our system is described by a quantum field theory and the 
existence of a classical configuration of domains requires a quantum-to-classical 
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Figure 3.2: Plot of the Fourier-space structure factor, S{k,t) = k'^G{k,t), vs. co- 
moving momentum k at various times during the evolution for quench parameter 
r = 1. 



transition, of which decoherence is an essential condition. In a decohered system 
we could, in principle, extract a positive-definite probability distribution of field 
configurations, -P[$], from the density matrix. Classical defect solutions would 
then appear as particular field configurations drawn from [0. 

In this work our focus is not on the extraction of classical defects from the 
quantum system but instead on the long wavelength modes that determine the 
size of correlated domains and, therefore, the average defect separation. The long 
wavelength modes interact with an environment of short wavelength and thermal 
fluctuations that destroys quantum coherence among the long wavelength modes 
and results in a finite correlation length for the system While quantum 

fluctuations are indeed important to the description of phenomena within the 
scale of one domain, beyond this scale we may hope to treat the modes as classical. 
With this abbreviated estimate of decoherence we show that the existence and 
size of correlated domains is indicated by an infrared peak in the power spectrum 
of the equal-time momentum-space two-point function. 

For a free field theory quenched into the unstable region by an instantaneous 
change in the sign of the square of the mass at t = to it is possible to obtain an 



analytic expression for the equal-time momentum space two-point function |^ 
After the quench, the momentum space structure function k'^G{k,t) has a strong 
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maximum at the value 



kmax-^^j"^. (3.45) 

The real space Fourier transform of this function (normalized to unity at |r — Fl = 
0) is 

Gdf-r'l.t) ~ e 8^ . (3.46) 

This form of the equal-time correlation function is common in condensed matter 
systems . From the equal-time real space correlation function we can identify 




the correlation length 



W— -7-^- (3.47) 

"'max \^) 

In analogy with the free field theory, we assume the domain size is proportional 
to the size of the maximum of the momentum space structure function. 

Another estimate of the domain size is motivated by the counting of defects 
in a classical condensed matter system. While the identification of topological 
defect structures can be a complicated task, when defects are well-formed (so 
that the width of the defect is much smaller than the typical defect spacing), it 
is possible to identify zeros of the classical field configurations with topological 
defects. In a model with a global Oin) symmetry in n spatial dimensions, and 
when the field probability distribution is Gaussian, a formula for the ensemble 
average density of field zeros was given by Halperin and derived explicitly by 



Liu and Mazenko 1167 



t^l — — , <-^2 — 7^, t-'S — ^ 



For a single scalar field in 3 spatial dimensions Eq. ( ^.48| ) with n = 1 is valid as the 
ensemble- averaged density of zeros along a one- dimensional section of the field. 
The validity of the Gaussian approximation is further discussed in A zero 
of the classical field configuration does not uniquely identify a topological defect. 
Thermal fluctuations give a large number of zeros of the field configuration on 
small scales |^|. The zeros of the field configuration which are not associated with 



different defects lead to an ultraviolet momentum divergence in the expression 
for the zero density, Eq. ( p.48| ). If we are to make physical sense of the zero 
formula, a coarse graining of the field configuration is needed. To effect this 
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coarse-graining in Eq. ( p.48| ) we impose a spatial momentum cut-off at tlie upper 



edge of the unstable momentum band. With this choice only the unstable modes 
will contribute to average density of defects, p{t). The length scale set by the 
coarse-grained density of defects provides another estimate of the average domain 
size. 

When the structure factor k'^G{k,t) is very strongly peaked about kmax the 
length scales set by kmaxif) and pit) differ only by a numerical factor of order 
unity. In this case we can approximate the unstable portion of the spectrum as 

k^G{k,t)r^6{k-k^a.)- (3.49) 
The zero density for a single scalar field in 3 dimensions is then 

1 f ik^Gik,tW'^ , , 

Figure |3.3| affirms this single-scale relation using the linear model introduced in 



Sec. p.4.2| . Additional evidence for a one-scale model was given in [|TOl in the 
context of the dynamics of a classical field phase transition in 1 + 1 dimensions. 
For very weak coupling not considered here and when a peak is no longer evident 
in S{k,t), the one-scale approximation will fail and it is possible in principle 
for the defect density as measured by p{t) to depart significantly from one per 
correlation volume [ |7T| . 



3.4 Critical scaling 

The evolution of the two-point function and the extraction of correlated domains 
allows a first-principles analysis of the mechanisms important for topological de- 
fect formation. In light of the freeze-out scenario it is interesting to compare the 
initial size of correlated domains for scale factors a^-it) with different expansion 
rates r. To quantify the dependence of the initial size of correlated domains on 
the quench rate r of the phase transition, we compare the domain size for dif- 
ferent values of the parameter r. As discussed above, the average domain size is 
proportional to the maximum in the infrared portion of S{k, t). We compare the 
maximum kmax for different values of r and at two distinct sets of times during 
the phase transition. The first set of domains is measured when the square of the 
effective mass reaches a minimum value. This provides an early-time measure of 
the size of domains. The results are shown in Fig. |3.4| . Filled squares are mea- 
surements from the numerical simulations. The solid line is a plot of the best-fit 
linear function to the data, 

logio(A;„.a.) = -0.351ogio(r) - 0.14. (3.51) 
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Figure 3.3: Plot of the ratio R{t) = vs. cosmological time t for quench 



parameter r = 10 using the early-time model of Sec. p. 4. 2 
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Figure 3.4: Log-Log plot of k^ax vs. r in the case of slow, underdamped quenches 
(r > 1.0) for domains formed early in the phase transition, at the time when the 
square of the effective mass reaches a minimum value. 
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Figure 3.5: Log-Log plot of kmax vs. r in the case of slow, underdamped quenches 
(r > 1.0) for domains formed late in the phase transition, at the time when the 
square of the effective mass reaches a local maximum value. 



The second set of domains is measured when the square of the effective mass 
reaches a local maximum value after the phase transition. This provides a late- 
time measure of the size of domains. The results are shown in Fig. ^]5| and 



Fig. 3.6. Filled squares are measurements from the numerical simulations while 



solid lines plot the best-fit linear function to the data. In Fig. |3.5| the best-fit 
power-law is 



logio(fcri 



while in Fig. 3.6 the best fit is 



logio(fc« 



-0.35 logio(r) -0.26. 



-0.28 logio(r) -0.22. 



(3.52) 



(3.53) 



For slow quench rates (r > 1.0), at both early and late times, the dependence 
of kmax on the quench rate r is well approximated by a power law. 



r 



-0.35 



(3.54) 



For fast quench rates (r < 1.0), the dependence of k^ax on the quench rate r is 
well approximated by a power law. 



kr. 



r ~ r 



-0.28 



(3.55) 



The origin of these exponents and the difference between slow and fast quenches 
is discussed in the following section. 
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Figure 3.6: Log-Log plot of kmax vs. r in the case of fast, overdamped quenches 
(r < 1.0) for domains formed late in the phase transition, at the time when the 
square of the effective mass reaches a local maximum value. 



3.4.1 The freeze-out scenario 

In the Zurek scenario the frozen correlation length and therefore the initial size 
of correlated domains scale as a power law of the quench time r as 

^freeze ~ (3.56) 

where /x and u are the equilibrium critical exponents for the correlation length 
and relaxation time respectively. The correlation length and relaxation time 
are identified as, respectively, the length and time scales that characterize the 
equilibrium behavior of the propagator near the critical temperature |7^. Under 
the scaling hypothesis, the equal-time propagator is written 

Gk = k-^^-''^F[h{e)k]. (3.57) 

Here e is the reduced temperature, Eq. (p^), is a critical exponent (not to be 
confused with conformal time) and F is a dimensionless function. The finite- 
temperature equilibrium correlation length is 

^(e) = h{e) ~ e-". (3.58) 

Similarly the two-time propagator is 

G^^k = cu°F[/i'(e)c<;, h{e)k] (3.59) 
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so that the relaxation time 



r h' e ^. 

In the Hartree-Fock approximation the equihbrium propagator is 



(3.60) 



1 



1^2 _ f^2 _ ml^f 



coth( 



eff 



(3.61) 



where the effective mass is given by the equihbrium value Eq. ( p.26|) . Near the 
critical point 



m 



(3.62) 



where A is a constant independent of temperature. The scaling behavior of the 
propagator is as and ^ and therefore the critical exponents for the theory are 



1 
2' 



H = V 

In the Zurek scenario we therefore expect the power law scaling, 



~ T 3 . 



(3.63) 



(3.64) 



Our observation of the power law exponent of 0.35 in the numerical simulations 
for slow quenches, r > 1.0, (Fig. 



and Fig. p.5|) is consistent with this result. 
The calculation of the equilibrium critical exponent for the relaxation time 
depends on the dynamics of the mode functions for the the low k modes near the 
critical point and in particular whether the dynamics is overdamped or under- 
damped. Near the critical point, ^'ijf ~ 0, and the mode function equation is 
approximately 



d? aitc) d k'^ 
dt'^ a{tc) dt a?{tc 



fk{t) = 0. 



(3.65) 



The critical time tc can be estimated as the time when the equilibrium effective 
mass, Eq. ( p.92| ), goes to zero. 



(3.66) 



and the value of the scale factor at tr is 



a{tc 



'tc + T 



T 



1.3. 



(3.67) 
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The mode function equation near the critical point Eq. ( p.65[ ) is that of a damped 
harmonic oscillator with natural frequency Uq = and damping constant T 
where 

r = 3H{t,) ^ -. (3.68) 
r 

The modes for which Uq > T or 

1-3 

k> — (3.69) 
r 

are underdamped. The range of wavenumbers k for the modes that determine 
the size of correlated domains also depends on r. Eq. (|3.69| ) is a condition on r 
such that when 



r > r, (3.70) 

the modes responsible for domain formation are underdamped and r^, is to be 
determined. For the underdamped case we expect fc^ax ~ t^^. Using k^ax 
as a representative wavevector of the modes responsible for domains, Eq. ( |3.69|) 
implies the underdamped condition 

n ~ 1. (3.71) 

Slow quenches with r > 1.0 are therefore underdamped. Fast quenches are 
overdamped and the dynamics of the mode functions will be dominated by the 
first time derivative. In the overdamped case, the critical exponent, /z, now 
assumes the non-relativistic value 



overdamped founder damped: (3.72) 

and 



kmax ~ r 4 (3.73) 



This is consistent with the power law exponent of 0.28 measured for overdamped, 
fast quenches (r < 1.0) shown in Fig. |3.6| . In both the overdamped and un- 
derdamped cases, the prediction of the "freeze-out" proposal for the power law 
exponent of the scaling of the initial defect density with the quench rate appears 
to be in excellent agreement with the numerical simulations. 



3.4.2 Early-Time linear approximation 

If the initial size of correlated domains and the initial density of topological 
defects are determined by processes occurring very near the critical point, as 
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Figure 3.7: Log-Log plot of kmax vs. r in the analytic model, Eq. ( p.82| ). The 
maximum kmax was determined as soon as a peak was evident in the structure 
function S{k, t). 



claimed in the freeze-out proposal, then the power law scaling of domains with 
the quench parameter r observed in the simulation of the Hartree-Fock mode 
function Eq. ( p.23|) will also appear in a linear approximation. A linear amplitude 
approximation to the full mode function equation is valid for slow quenches and 
only near the onset of the spinodal instability, before the unstable modes have 
grown appreciably and back reaction is important. 
We consider the linear equation 



^ + 3-- + ^ 



dt 



a dt 



m 



0, 



(3.74) 



where 



K' = e + ki 



(3.75) 



and fco is related to the initial temperature. 



''o- 24 °- 



(3.76) 



In conformal time we have 



4r (t + r) 



(3.77) 
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with conformal modes fk defined by 



fk = Fktt. 

Equation ( |3.74} ) may be rewritten in conformal time as 



fk = 0. 



(3.78) 



(3.79) 



Tfiis equation can be solved exactly in terms of parabolic cylinder functions. 
However, a simpler solution is obtained for times near the critical point. We 
further approximate by expanding around the critical point, 



d"^ mK , ' 



fk = 0, 
2Kt 



m 



(3.80) 
(3.81) 



The properly normalized solution with vacuum boundary conditions before the 
instability is 



fk 




2 / mKx^ 
3 



r 



where 



X 



(3.82) 



(3.83) 



and Ki/s is a modified Bessel function. We now use the analytic solution Eq. ( p.82|) 



to examine the dependence of the position of the peak in the Fourier space struc- 
ture factor k'^G{k, t) on the quench parameter r. The location of this peak kmax{t) 
redshifts throughout the phase transition, moving towards lower momentum as 
the domains coarsen. In order to compare domains formed in the linear model 
with those of the Hartree evolution we compare the position of the peaks at very 
early times when a peak is first identifiable in S{k,t). Plots of S{k,t) with the 



analytic modes were made with Mathematica. The results are shown in Fig. pTZ . 
Filled squares are measurements from plots of the analytic modes. The solid line 
is a plot of the best fit linear function to the data. 



logio(fc« 



-0.34 logio(r) -0.09. 



(3.84) 



The scaling of domains with the quench parameter r with the same power law 
exponent as seen in the full numerical simulation for slow quenches is evident. 
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3.5 Summary and future work 



Using a two-loop truncation of the 2PI-CTP equations of motion for the two- 
point function of a quantum scalar field undergoing a phase transition in a 3-1-1 
dimensional spatially flat, radiation-dominated FRW universe we have shown 
that the size of correlated domains, measured as the maximum of the peak of the 
infrared part of the spectrum of k'^G{k, t), scales as a power of the quench rate r 
as 

idamains ~ T°'^^, (3.85) 

for slow, underdamped quenches (r > 1.0) and 

idcmiains ^ T°'^^, (3.86) 

for fast, overdamped quenches (r < 1.0). The observed power-law scaling of 
correlated domains is quantitatively consistent with the freeze-out hypothesis. In 
both overdamped and underdamped cases, the value of the power-law exponent 
extracted from evolution of the two-point function is in good agreement with the 
value calculated in the freeze-out scenario using the critical scaling exponents for 
a theory in the Hartree-Fock approximation. 

To further explore the behavior of the quantum system near the critical point 
wc introduced an approximate linear model valid for slow quenches and short 
times after the onset of spinodal instability. In this linear model, the size of 
correlated domains scales with the quench rate r to the same power as observed in 
the underdamped numerical simulations and predicted by the freeze-out proposal. 
This provides analytical evidence for the freeze-out hypothesis. 

In order to place these results in physical context it is important to understand 
the validity of the two-loop truncation of the 2PI-CTP effective action. The 
truncation of the effective action to any finite loop order is an approximation. On 
some time scale we expect the contribution from higher loop terms to modify the 
two-loop dynamics. Unfortunately, the higher-loop equations are time-nonlocal 
and an exact determination of this time scale is not possible. Wc provide instead 
an approximate analysis. The higher-loop terms in the 2PI-CTP effective action 
contain additional powers of the coupling constant A. In the perturbative regime 
and in the small-coupling limit these terms are suppressed. This is the case early 
in the evolution, when the effective mass is positive. However, the dynamics 
of the second-order phase transition is dominated by spinodal instabilities, the 
exponential growth of long-wavelength fluctuations visible as a peak in k'^G{k, t). 
At times slightly beyond the onset of the phase transition, the effective mass 
is negative and long wavelength modes are growing. At these early times the 
growing modes have yet to affect the dynamics and higher-loop terms remain 
negligible. Eventually, the growth of fluctuations causes the effective mass to 
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grow until they dominate at the spinodal point (when the effective mass reaches 
zero from below), 

A J dkPG{k,t) r^m^. (3.87) 

Since, for example, the three-loop piece of the effective action is 

Ts-zoop ~ A^G^, (3.88) 

higher loop terms are likely to be important beyond the spinodal point. Whether 
these terms actually alter the size of domains during the late stages of the phase 
transition is an open question. The two-loop truncation of the 2PI-CTP effec- 
tive action, in the absence of a mean field, is also exact to leading order of a 
approximation. Although there are no defects in the leading-order Large-N 
model, domains can form in a symmetry-breaking phase transition. While the 
interpretation of our model as a large-N approximation provides mathematical 
rigor, it is not clear to what physical system the large-N approximation applies. 

In spite of these difficulties our approach remains, to date, the most com- 
prehensive treatment of the dynamical formation of domains in nonequilibrium 
quantum field theory. In contrast to previous work we have allowed for the process 
of back reaction, which permits the quantum field to exit the region of spinodal 
instability. The power law exponent is the same whether it is measured at very 
early times with the analytic linear model (Fig. p.7|) or, in the numerical simula- 
tions, at early times in the middle of the spinodal region at the minimum value 
of TTT'lff (Fig- P-4| ) or, at late times in the stable region at the maximum value of 
m^jj (Fig. p.5|) . This is consistent with the idea that the relevant processes for the 
growth of domains occur near the critical point. However, back reaction is essen- 
tial to the freezing in of the value of the power-law exponent and to the freeze-out 
hypothesis. Although the power-law is accurately recorded in the early-time an- 
alytical model, if back reaction is ignored and the linear model is (incorrectly) 
extrapolated to late times, the predicted scaling exponent is different 

A;max~r-i (3.89) 

The results derived from the numerical simulation of the mode function equa- 
tion with back reaction and the analytical solution around the critical point are 
evidence supporting the first verification of the freeze-out scenario in 3 spatial 
dimensions and in a realistic system relevant to the early universe. It is at first 
surprising that arguments based on classical equilibrium critical scaling apply in 
the context of dynamical quantum field theory. A possible explanation rests with 
the high temperature initial conditions and the properties of classical A$^ theory 
in FRW spacetime. Since the initial temperature is much higher than the initial 
effective mass 

To » me//(0) (3.90) 



45 



0.0 2.0 4.0 6.0 8.0 10.0 12.0 14.0 16.0 
t 



Figure 3.8: Plot of the square of the effective mass M = fri^ff vs. cosmological 
time t in the equihbrium case (dashed hue) where "^e// ~ —m^ + 24^?ft) ^^'^ 
in nonequihbrium case (sohd hne) where '^e// determined by the full numerical 
simulations. 



the field is approximately conformally invariant. In an FRW universe, a confor- 
mally invariant field initially in thermal equilibrium will remain in equilibrium 
at a redshifted temperature 

m = (3.91) 

A well-known example is the redshifting of the blackbody spectrum of the cos- 
mic microwave background radiation as the universe expands. To illustrate the 
approximate conformal invariance in our model we compare the dynamical ef- 
fective mass from the simulations with the effective mass of a theory with the 
equilibrium redshifted temperature of Eq. (|3.91|) 



The results are shown in Fig. |3.8| . The slight difference in the equilibrium and 
nonequihbrium curves at the beginning of the evolution is due to the difference 
between the general Hartree-Fock effective mass given by Eq. (|3.25| ) and the 
high temperature, small A limit given by Eq. ( p.26| ). Only after the onset of 
the spinodal instability do the equilibrium and dynamical effective mass differ 
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significantly. The approximate conformal invariance of tlie tlieory means tliat 
tlie finite-temperature effective potential and the critical behavior derived from 
it are approximately valid until times near the onset of the spinodal instability. 
The high temperature initial conditions also offer a possible explanation for the 
observed classical behavior of the system. At such high temperature, thermal 
fluctuations dominate over quantum vacuum fluctuations. 

The agreement in the value of the power law exponent between the micro- 
scopic evolution equations of the quantum field theory and phenomenological 
critical scaling supports the contention that quantum critical systems in the early 
universe share, in certain circumstances, many of the properties of their classical 
counterparts. It would be very interesting to further explore these connections. 
To do so, however, it is necessary to go beyond the (relatively) strong coupling 
and high temperature conditions used in this research. 

The results derived here represent a preliminary step in a first-principles ap- 
proach to the calculation of the topological defect density immediately following 
the completion of a second-order phase transition in the early Universe. Our 
work uses a microscopic quantum field theory and incorporates realistic initial 
conditions and a physical quench mechanism. However, the domain wall defects 
formed in this model are inconsistent with cosmological observations. It is not 
difficult in principle to apply the methods used in this paper to more cosmologi- 
cally realistic theories such as cosmic string models. The analytical and numerical 
evidence for the power-law scaling of the domains is expected to hold in a more 
realistic model. 

The formalism of the 2PI-CTP effective action and related techniques can be 
used to probe more general questions related to the physical aspects of quantum 
critical dynamics |7^. To observe and isolate quantum processes it is necessary 
to relax the assumption of high temperature initial conditions, allowing quantum 
vacuum fluctuations to dominate over thermal fluctuations. Also important are 
the detailed properties of the system-bath interaction. A system-bath interaction 
such as proposed in could address such issues as critical slowing down and 
the role of dissipation and noise. Work is in progress along these lines. 
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Chapter 4 



Classical critical dynamics: Formation and 
evolution of topological textures 



4.1 Introduction 

In Chapter 3 we used the techniques of nonequihbrium quantum field theory to 
examine the dynamics of a phase transition in the early universe. Unfortunately, 
the approximation methods necessary to obtain even a numerically soluble quan- 
tum dynamical system preclude a detailed look at topological defects themselves. 
In this chapter, that in part represents work published in [^, we examine the 



phase transitions of a classical field theory in which topological defects emerge 
as localized, stable solutions. Since the full dynamics of a classical field can be 
simulated numerically, this approach offers another window into the nonperturba- 
tive processes important for the formation and interaction of topological defects. 
Specifically, we study the formation of topological textures in a classical 0(3) 
scalar field theory in 2 + 1 dimensions. Systems containing topological textures 
are of particular interest as they provide an arena in which to explore myriad dy- 
namical issues ranging from nonequihbrium phase transitions to the coarsening 
dynamics of quantum field theory. 

Topological textures form in d spatial dimensions when the vacuum manifold 
M is such that nrf(M) is nontrivial. Unlike other topological defects, textured 
fields never leave the vacuum manifold but instead are "knots" of gradient en- 
ergy, given relative stability by their topological winding around the vacuum. 
Like cosmic strings, global textures have been studied extensively as candidates 
for generators of large-scale structure formation |7^, In condensed matter, 
topological textures are studied in a wide variety of (mostly two-dimensional) sys- 
tems, from superconductors |76| and quantum hall ferromagnets to superfiuid 
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Nonequilibrium phase transitions and topological textures 



The classical dynamics of textured systems provides an ideal study for nonequi- 
librium phase transitions. The richness of texture interactions involves multiple 
dynamical length scales. Unlike other topological defects, textures do not have a 
fixed size and the evolution of a textured system is enriched both by the interac- 
tions between textures and by the changing scale of the texture itself. In 2 + 1 
dimensions it is observed that the phase-ordering of a system containing topo- 
logical textures is characterized at late times by at least three length scales: the 
(average) texture size, texture-texture separation and texture-antitexture sepa- 
ration 1^9], In distinction, the coarsening dynamics of systems containing 



domain walls, vortices or monopoles are effectively single-scale and described, at 
late-times, only by the average defect separation. In a nonequilibrium phase tran- 
sition the multiple independent length scales of textured systems are determined 
by the phase transition dynamics. However, the arguments of the Kibble-Zurek 
mechanism are ambiguous when applied to textured systems. Which, if any, of 
the multiple length scales observed in a texture distribution should we equate 
with the single frozen correlation length? Although the dynamics of a textured 
system can be qualitatively different, the details of the formation of topological 
textures in a nonequilibrium phase transition are mostly unexplored (see however 
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Topological textures and the large-N expansion 

A detailed understanding of the field dynamics of a classical textured system 
can also be used to better understand the large-N expansion, a common tool in 
the study of nonequilibrium quantum field theory (see e.g. and references 
therein). The large-N expansion truncates both the dynamic and quantum con- 
tent of the full quantum field theory by replacing the physical system under study 
with one containing a very large number of scaler fields. In doing so, the dynam- 
ics of the quantum system are dominated by the dynamics of Goldstone modes 
and the density matrix is truncated to Gaussian form [^. However, real quan- 
tum systems contain only a small number of fields and are described by density 
matrices that may be far from Gaussian. It is therefore important to understand 
to what extent the large-N approximation reflects the physical processes of a real 
system. 

Phase ordering in d spatial dimensions of > c? + 1 scalar fields provides a 
particular example of a useful theoretical laboratory in which to analyze the dy- 
namical and quantum content of the large-N approximation method. If A^ = d+1, 
the dynamics of phase ordering is the dynamics of interacting topological textures. 
If N > d + 1, no topological defects are formed. In either case, phase ordering 
is driven by massless Goldstone dynamics. However, the Goldstone dynamics 
of large-N can be radically different from a real system. As we discuss in the 



49 



next section, the coarsening of a 2 + 1 dimensional system containing topolog- 
ical textures is observed to depend on three separate length scales which grow 
with different powers of time. A large-N analysis of the same system misses the 
contribution of topological textures and predicts only one simple scaling length 

J^large-N ~ ^2 . 

To understand how the large-N approximation fails, the full classical equations 
of motion can be studied numerically. For example, assume the dynamics of the 
coarsening system is described by an N-dimensional order parameter $ with fully 
dissipative time-dependent Landau- Ginzburg equation (TDGL), 

ai>(x,t) _ 5Fm 

dt 5<i> ' ^ ^ 

with free energy 

= j - rrv'f]. (4.2) 

The dynamical content of the large-N approximation is revealed as this system is 
solved for larger and larger values of N. In the formal limit of infinite N, both the 
classical and quantum systems can be solved exactly (albeit numerically). In the 
infinite-N case, the comparison of classical and quantum ordering dynamics yields 
insight into the true quantum content of the large-N approximation. The focus in 
this current work is on the formation of topological textures in a nonequilibrium 
phase transition of a classical field theory. We hope to return to the role of 
textures in the large-N approximation and nonequilibrium quantum field theory 
in future reports. 



Organization 



This chapter is organized as follows. Section reviews the properties of topo- 
logical textures in 2 -|- 1 dimensions and the scaling violations observed in textured 
systems when quenched from a disordered phase. Section describes two mod- 
els of a nonequilibrium phase transition, a dissipative quench and an external 
quench. In addition we provide the numerical techniques used in the evolution 
and the length scales used to characterize the texture distribution. Section |4.4| 
contains the main results of this chapter: an explanation of the quench rate 
dependence of the average texture separation and the average texture width ob- 
served near the end of the external quench phase transition. We show that the 
Kibble-Zurek mechanism is recovered at early times but that by the end of the 
phase transition, formation dynamics and phase ordering dynamics are intrinsi- 
cally linked. We use these results to argue that textured systems carry an imprint 
of the freeze-out correlation length to late times and suggest the possibility of 
novel, late-time measurements of dynamical critical phenomena. Concluding re- 



marks are made in Sec. 4.5. 
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4.2 Topological textures in two spatial dimen- 
sions 



The properties of systems containing topological textures differ in many ways 
from those containing other topological defects. To illustrate these differences, 
consider an 0(3) invariant model with classical action, 

S = J d^x[d^$ ■ + ^(|$p - m^], (4.3) 

and resulting equations of motion 

- v'$. + A$.($' - m^) = 0. (4.4) 

For the classical dynamics that is of interest here, the value of the quartic self- 
coupling A and the mass parameter may be scaled to unity with an appropriate 
scaling of the field and spacetime units, 

$4 ^i'm'^, t — i> t^/Xm, X — > x\f\m. (4-5) 

The vacuum manifold is characterized by $ ■ $ = 1, topologically a 2-sphere. 
When the field is constrained to the vacuum, the action is that of a nonlinear 
sigma model (NLSM). For asymptotically uniform fields, a vacuum field configu- 
ration is a map 

■ real space fieldspacei \ / 

divided into distinct topological sectors by winding number n ||8^ . In two spatial 
dimensions these windings are topological textures, topologically stable and static 
solutions with finite energy. For example, the following field configuration, a 
texture of negative unit winding, wraps around the vacuum 2-sphere exactly 
once, 

Aarsmid)) ^ Aarcosid)) ^ r"^ — Aa^ , , 

$1 = — $0 = — = (A 7) 

1 r3 + 4a2 ' ^ r2 + 4a2 ' ^ + Aa^' ^ ' ^ 

At the origin, $ points in the —$3 direction while for r — > cxo it points in the 
opposite direction, $3. In between, at r = 2a, the field points radially outward 
like a hedgehog. The parameter a characterizes the size of the texture and, 
in contrast to other defects such as strings or monopoles, is independent of the 
parameters of the action Eq. ( [4.3|) . The NLSM also possesses an exactly conserved 
topological charge Q, which is expressed as the integral over a topological charge 
density p, 

p = ■ (dj A dy$) (4.8) 
47r 



Q= / d'xpix,y) (4.9) 
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For the single winding texture, Eq. ([4.7[), the topological charge density has a 
particularly simple form 

1 4a2 , , 

P = , ^ , , 4.10 

with total topological charge Q = —1. The energy E = ^ J d?x{d^i ■ d^i) of 
the single winding configuration Eq. ( [4.7|) is simply E = Sir, independent of the 
size a of the texture. In accordance with Derricks theorem the energy of 
a texture configuration in higher dimensions scales as a positive power of its 
size, demonstrating their instability to collapse. In lower dimensions textures are 
unstable to expansion. The stability of a texture configuration in any dimension 
can be altered by adding higher derivative terms to action, as was done in the 



Skyrme model of nucleons p5 



Although the energy of a single isolated topological texture is independent of 
its size, systems with multiple textures are not static. Textured systems order 
under the constraint of conserved topological charge Q. However, the details of 
texture interactions are not well understood. It is observed that textures and 
antitextures can annihilate with each other and that more isolated textures can 
decay by unwinding |]7D|, ^ . The details of the unwinding process depend on the 



model under consideration. For the action Eq. (13) textures can decay by pulling 
the field off the vacuum manifold and unwinding. This is analogous to isolated 
textures in higher dimensions which are unstable to shrinking and unwind when 
their gradient energy is large enough to move the fields off the vacuum manifold. 
In a NLSM where the fields must remain on the vacuum manifold, textures can 
unwind through higher-derivative terms present in the lattice discretization. 

Topological textures and the violation of dynamical scaling 

When a disordered system is quenched into the ordered phase (e.g. by the removal 
of thermal fluctuations), the approach to a new equilibrium and associated long- 
range order is through a sequence of nonequilibrium states. A simple example is 
the Ising ferromagnet in 3 spatial dimension. Above the Curie temperature, spins 
are randomly oriented and there is no net magnetization. A snapshot immea- 
diately following the quench of this system to zero temperature shows a system 
containing many domains in which the spins are coherently aligned. The net 
magnetization remains zero because the domains are randomly oriented with re- 
spect to each other. However, the true T = ground state of this system consists 
of only one domain and no domain walls. The system approaches this equilibrium 
through the dynamic process of coarsening, during which small domains shrink 
and large domains grow. The process of coarsening has been studied extensively 
in both condensed matter systems and in the early universe (for reviews see |Q 



and |36|). The scaling hypothesis arose from these studies and states that, at 
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late times, and in a quench from a disordered state, the domain distribution is 
statistically characterized by a single dynamical scale L{t) which grows in time. 
For simple dissipative systems, the approach to (T = 0) equihbrium is governed 
by model-A time-dependent Landau- Ginzburg (TDGL) dynamics for a noncon- 
served order parameter 

dt 5$ ■ ^ ' ' 

It has been shown both theoretically and experimentally that the average domain 
size L{t) grows with time as 

m = (4.12) 

which is also the power-law predicted by a dimensional analysis of Eq. ( [4. 111 ). 
Not all systems coarsen in such a simple way and there is currently no general 
theoretical framework in which to explain why some systems scale while others 
do not. In particular, systems containing topological textures in one and two 



spatial dimensions are strong exceptions to the scaling hypothesis fTy, ^ In 



a 2 + 1 dimensional 0(3) model with dissipative dynamics and an instantaneous 
quench from the disordered phase, late-time coarsening of topological textures is 
described by at least three different length scales. These three length scales are 
the average defect-defect separation L^ep, the average defect width L^, and the 
average texture- antitexture separation Ltat- The scaling hypothesis is violated 
since at late times these three length scales are observed to grow with different 
powers of time 

Lsep{t) = ^ItK (4.13) 

LUt) = (4.14) 
Ltat{t) = (4.15) 

where is the correlation length in the disordered phase. It is remarkable that 
length scales describing the late-time ordering of textures retain a dependence 
on the initial correlation length. In distinction, systems that obey the scaling 
hypothesis erase any memory of their initial state and at late times retain only a 
dynamical scale (for example Eq. ([4.12|) in the case of an Ising ferromagnet). The 
scaling of systems that contain topological defects such as vortices implies that 
the vortex distribution at late times is independent of the dynamical details of 
the phase transition. For textured systems, it is not possible to have such a clean 
separation between the dynamics of the phase transition and the late-time texture 
distribution. In the context of a nonequilibrium phase transition we argue that 
is determined not from initial conditions but from the Kibble-Zurek mechanism, 
thus connecting early-time critical dynamics and late-time coarsening. In Sec. 
we discus how this connection may provide late-time probes of nonequilibrium 
critical dynamics. 
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4.3 Phase transitions and nonequilibrium clas- 
sical dynamics 

Classical interacting field theories appear as useful physical descriptions in diverse 
systems ranging from hydrodynamics and condensed matter to general relativity 
and the early universe. Since the dynamics of a classical field can often be sim- 
ulated numerically, classical approaches offer a window into the nonperturbative 
dynamics of the decohered sector of quantum field theory. In this section we de- 
scribe two models of the classical dynamics of a nonequilibrium phase transition, 
a dissipative quench and an external quench. It is useful, however, to first discuss 
the nature of a classical field in thermal equilibrium. 

A field contains an infinite number of degrees of freedom. For a field confined 
in a cavity, these excitations are the (countably) infinite number of standing 
waves, satisfying the cavity boundary conditions. If the field is maintained in 
equilibrium (for example by keeping the walls at constant temperature), each 
degree of freedom holds a finite average energy {E ^ T for a free field) given by the 
equipartition theorem. A classical field in a finite cavity in thermal equilibrium 
thus contains an infinite amount of energy! The energy density is also divergent. 
For black body electromagnetic radiation in a 3-dimensional cavity the spectral 
energy density Pt{^) is given by the Ray leigh- Jeans formula, 

pT{y)di' ~ Ti/dv. (4.16) 

The divergence of the spectral energy density at high frequency is known as 
the ultraviolet catastrophe and its resolution led Planck to the foundations of 
quantum theory. 

In light of the ultraviolet catastrophe, it what sense does equilibrium exist 
for a classical field theory? If the classical theory is only a long-wavelength ap- 
proximation to a quantum field theory, then a correct quantum accounting of the 
short-wavelength modes will suffice to render equilibrium meaningful. However, 
the classical Landau- Ginzburg field theories of interest here are a only contin- 
uum approximation to a fundamentally discrete condensed matter system. This 
discrete nature provides a natural cutoff in the frequencies of field oscillations, 
rendering the spectral energy density and total energy finite. In principle, the 
observables computed from a thermal classical field theory with a cutoff depend 
upon the value of cutoff. In practice, for certain situations and especially in the 
vicinity of a second-order phase transition, low frequency modes are remarkably 
independent of the behavior of the high frequency modes, and therefore insen- 
sitive to the cutoff. In this case, observables computed from long-wavelength 
modes are universal. 
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Figure 4.1: Plot of the topological charge density for a small section of the lattice 
at a time near the onset of the quench, t = 0.0. 



4.3.1 Dissipative quench 

As in Chapter 3, to study the dynamics of a phase transition, we are faced with the 
challenge of implementing a realistic quench. One possibility is to allow energy 
dissipation of the field to occur solely through a phenomenological dissipative 
term 7](j) in the equations of motion. The field evolves with equations of motion 

v''^'i + *^('5'-l) + #i = (4.17) 

and initial conditions are chosen so as to correspond to a high-temperature sym- 
metry restored phase. The initial conditions are implemented either through a 
Langevin equation or simply by assigning random values to (p and The ad- 
vantage of this approach is that the quench proceeds naturally. Energy is lost 
through dissipation but the dynamics of the fields are otherwise unaffected. In 
particular, and in distinction to other approaches in the literature |T^, no a 
priori form of the quench rate or the effective mass is assumed. 

Using the numerical techniques explained in Sec. |4.3.2|, details of the dynamics 



of the dissipative quench phase transition and the distribution of topological 
textures were extracted from a numerical solution of Eq. ([4.17|) . A series of 
snapshots of the topological charge density for the above dissipative quench model 
with rj = 0.8 is shown in Figs. |4.1| - |4.6| . Horizontal axes label location on the spatial 



lattice while the vertical axis labels the value of the topological charge density. 
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Figure 4.2: Plot of the topological charge density at time t — 9.0. 




Figure 4.3: Plot of the topological charge density at time t — 10.5. 
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Figure 4.4: Plot of the topological charge density at time t — 12.0. 




Figure 4.5: Plot of the topological charge density at time t — 13.5. 
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Figure 4.6: Plot of the topological charge density at time t = 17.5. 



Though the 2 + 1 dimensional theory defined by Eq. ( |4.17| ) has no finite- 
temperature equilibrium phase transition, there can still be an effective dy- 
namical transition. As energy is dissipated, the high temperature initial state 
dominated by short wavelength thermal fiuctuations eventually evolves into a 
state dominated by long wavelength topological textures (well-defined peaks in 
the topological charge density). Textures become quasi-stable excitations when 
the effective temperature drops substantially below the texture energy of Svr. 
The dynamics of the phase transition is also evident in the time evolution of 
the total number of textures Q and of the order parameter, = (|$|), where 
brackets denote an average over the lattice. In Figs. and we exhibit the 
total number of textures and the order parameter for a quench with rj = 0.7. In 
the early stages of the phase transition both Q and decrease exponentially as 
dissipation dampens the fields. However, when enough energy is removed, the 
potential term in Eq. ( ^.17| ) becomes important and the system evolves to the 
stable ground state, = 1.0. Texture formation occurs primarily in the interval 
from t = 9.0 to t = 12.0. 

The dissipative quench model has one free parameter, the dissipation constant 
rj. In Fig. O we show the number of textures measured at the local maximum 
Qmax for different values of rj. Larger r] results in a faster quench. Surprisingly, 
it appears that, contrary to the Kibble-Zurek mechanism, faster quenches lead to 
a smaller number of topological defects. Unfortunately, we have yet to develop a 
convincing explanation for this behavior. It is possible that since this system has 
no equilibrium phase structure the arguments of the Kibble-Zurek mechanism 
simply do not apply. It is also possible that the identification of textures through 
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Figure 4.7: Plot of the total number of textures Q vs. time t for a quench with 
T] = 0.7 
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Figure 4.8: Plot of the order parameter (|$|)iat vs. time t for a quench with 
ri = 0.7 
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Figure 4.9: Plot of the total number of textures Qmax vs. dissipation constant t] 



the topological charge density is inadequate: we may mistake thermal fluctuations 
for stable defects. It was our original intent to study the Kibble-Zurek mechanism 
in the context of the multiple length scales of a textured system. In the following 
section, motivated by similar models in the literature [|l^, we introduce and 
analyze a 2 + 1 dimensional texture model with a phase transition implemented 
through an imposed equilibrium phase structure. The dissipative quench provides 
a useful physical model to which we hope to return our attention in future work. 



4.3.2 External quench 

In order to study the dynamical formation of topological textures in the absence 
of strong thermal fluctuations, the equations of motion Eq. ( |4.4|) are modified 
by the addition of an externally controlled time-dependent effective mass and a 
stochastic environment. The effective mass is 

mlff{t) = tanh (^^^^ , (4-18) 

where r is the quench rate and defines the critical point where spinodal insta- 
bilities, characteristic of a second order phase transition, begin to grow. Both in 
the early universe and in the laboratory, texture formation occurs in an environ- 
ment containing many other degrees of freedom. The effect of the environment is 
assumed to be that of a stochastic (Gaussian white noise) driving term ^{x, t) and 
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simple ohmic dissipation, obeying a fluctuation-dissipation relation. The external 
quench equations of motion are 

$i - v'-^-^ - + - m%f) = t), (4.19) 

where the fluctuation-dissipation relation is 

(e(x, t)i{x, t')) = 2T7]6\x - x)6{t - t'). (4.20) 

The fluctuations, controlled by the parameter T are weak and serve only to give 
a random kick as the flelds roll to the true vacuum. The dissipative term ?7$j 
is strong for short wavelengths and helps desensitize the classical fleld model to 
modes near the lattice cutoff (used below in the numerical simulations). Indepen- 
dence from the cutoff was also verifled operationally by halving the lattice spacing 
and observing no change in the output. The use of a time-dependent mass in the 
equations of motion of this 2-1-1 dimensional model deserves comment. At t = 
the effective mass is positive and the system fluctuates around the false vacuum 
|$| = 0. When t reaches tc the system is destabilized and the flelds begin to roll 
to the true vacuum |$| = 1. This behavior mimics that of a phase transition in 
3-1-1 spacetime dimensions. It is well-known that there is no continuous symme- 
try breaking and accompanying long range order in spatial dimensions of two or 
less . Therefore the time dependence of the effective mass does not come from 



changing a thermodynamic variable, such as temperature. Instead, we use the 
external time-dependent effective mass to allow for a controlled passage between 
the fluctuation-dominated and texture-dominated phases. Within this controlled 
setting the role of the quench rate in the formation of texture can be isolated and 
studied effectively. It is our hope that the lessons learned are applicable to more 
realistic phase transitions, both in the early universe and in the laboratory. 



Numerical parameters and techniques 



The dynamics of the non-linear stochastic system were solved numerically. Equa- 
tion ([4.19|) was discretized on a 2-dimensional lattice using a standard second- 
order leapfrog method for the time evolution |7D| and a second-order 
spatial discretization for the laplacian 



nn+ib]M 



l^n„b1W + ^(^(4b- + i]W-4Cb]W 



+$„[j-l][A;]+$„[j][A;-f l]-f $„b1[A: 
+mlff[ndt\)+U3\W\) 



l])-^n\j][ml[j\[k] 



(4.21) 



where 11 = ^ and x = V^t- The labels [j][k] identify the spatial lattice point 
and n labels the time step. At each time step and each lattice point the noise is 
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generated through a sum of M randomly distributed numbers 



161 



i=l ^ ' 

where —0.5 < < 0.5 are a set of A^^ random numbers. In the hmit 

when M — i> oo the distribution approaches that of a Gaussian with the 

variance required for the fluctuation-dissipation relation Eq. ( 4.20| ). In practice 



M = 24 was used and no change in the dynamics was observed for larger M. 
The timestep was dt = 0.02 and the lattice spacing was fixed at c/x = 0.5 on a 
square lattice with = 1000 sites per side. There are three physical parameters 
in the system, the quench rate r, the dissipation constant rj, and the strength 
of the fluctuations T. The quench rate was varied between r = 10 and r = 
150. Faster quenches (smaller r) deviate from the observed power-law behavior 
and asymptote to the behavior characteristic of an instantaneous quench as has 



been previously observed [|T8[. Longer quenches are in principle possible but 
take substantially more computer time. The dissipation constant was chosen at 
7] = 1.0 to produce relatively overdamped behavior. The evolution equations 
were supplemented with initial conditions generated by solving the equilibrium 
dynamics of Eq. ([4. 19]) with a value of the effective mass fixed at one quench 



timescale from the critical point, 

mljf{t < 0) = tanh (1.0). (4.23) 

Run times for the system to reach an order parameter of = 0.95 on a DEC 500 
Mhz workstation varied from hours for r = 10 to days for r = 150. 



Properties of the texture distribution 

Field configurations obtained from numerical simulations provide an embarrass- 
ment of riches. In order to fully explore the participation of topological textures 
in the dynamics of the phase transition, it is necessary to cull, from the abundance 
of information contained in the snapshots of the topological field configurations, 
a few well-chosen functions, whose time evolution provides a clear indication of 
the dynamics. The late-time coarsening of a textured system is described by at 
least three different length scales, the average defect-defect separation Lg^p, the 
average defect width L^, and the average texture- antitexture separation L^f. In 
the formation dynamics studied here we were able to measure reliably only Lgep 
and Luj. These length scales are defined through the the topological charge den- 
sity and the topological defect two-point correlation function. We consider the 
two-point defect density correlation function 

Ci(r,t) = (Ip(f)llp(f + f)|) - (|p(f)|)^ (4.24) 
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where () denotes average over the lattice. We define the normahzed correlation 
function 



(4.25) 

G{r, t) is a statistical measure of the "lumps" in the defect density function |p(a;)|. 
The average texture width Ly^ is proportional to the half-height scale r i of G{r, t) 
defined by 

G(r.,t) = l/2. (4.26) 
To measure the average texture separation we consider 

Qtotai = [ d^x\p\. (4.27) 



Qtotai counts the total number of topological defects in the system without dis- 
tinguishing between textures and antitextures. On dimensional grounds, 

Qtotai ~ -p— -^sj/s) (4-28) 

Sep 

where Lg^p measures the average separation of the topological defects and Lsys is 
the physical size of the lattice which did not vary between runs. These two length 
scales, and Lg^p, provide a detailed characterization of the topological texture 
network at the end of the transition. In principle, each scale offers a different 
window into the nonequilibrium dynamics of the phase transition. 



4.4 Unraveling critical dynamics 

In this section we explore the details of the external quench contained within 
the dynamics of Eq. ([4.19|) . Plots of the time evolution of the order parameter 
= ($^) and the average number of topological defects Q = f (Px \ p{x, t) \ are 
shown in Figs. [4.10| and [4.1 1| respectively. The plots correspond to an evolution 



with quench parameter r = 10. In the symmetric state {t < 10.0) both and 
Qtotai are close to zero. When t > 10.0 the external effective mass is negative 
and the fields, kicked by the small stochastic force ^{x,t), begin to roll to the 
true vacuum = 1. The order parameter is a measure of the number of 
lattice points where the fields have fallen to the true vacuum. As increases 
both the gradient energy and the number of defects increase. However as ^ 1 
most of the lattice is near the ground state vacuum and the defect formation 
process is over. Although the system is overdamped near freeze-out, the small 
oscillations in and Q near the end of the transition indicate that oscillations 
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Figure 4.10: Plot of the order parameter = versus time t for quench 

parameter t — 10. 
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Figure 4.11: Plot of the average number of textures Q{t) — J (Px\p{x,t)\ versus 
time t for quench parameter r = 10. 
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Figure 4.12: Log-Log plot of the average number of textures Q = J (fx\p{x,t)\ 
versus quench parameter r for different values of the order parameter = ($^). 
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Figure 4.13: Log-Log plot of the average texture width versus quench parameter 
T for different values of the order parameter (p = ($^). 
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around the true vacuum are not completely overdamped. After formation, the 
defect density does not remain constant but decreases slowly as texture-texture 
and texture-antitexture interactions dominate. 

Figs. and |]T| plot the average number of textures and the average 



texture width versus the quench rate. In each figure, symbols denote lattice 



measurements while solid lines denote best-fit power-laws. In Fig. [4.12| lattice 
measurements were made at early times (0 = 0.1) denoted by solid triangles on 
the lower graph, at intermediate times (0 = 0.5) denoted by stars on the middle 
graph and at late times (0 = 0.9) denoted by solid squares on the upper graph. 
The best-fit power-law for times early in the phase transition is 

log{Q)eariy = (1.31 ± 0.02) - (0.43 ± 0.01) log(r). (4.29) 
At late times the best-fit power law is 

log(Q)/ate = (3.17 ± 0.03) - (0.77 ± 0.02) log(r). (4.30) 
A crossover is apparent at intermediate times when the number of textures Q{t) 



is not given by a single power-law of the quench rate. In Fig. [4.13| lattice measure- 
ments were made at early times (0 = 0.1) denoted by solid triangles on the lower 
graph, at intermediate times (0 = 0.4) denoted by open circles on the middle 
graph and at late times (0 = 0.9) denoted by solid squares on the upper graph. 
The best-fit power-law for times early in the phase transition is 

logiLjeariy = (0.42 ± 0.02) + (0.14 ± 0.01) log(r). (4.31) 

At late times the best-fit power-law is 

logiLjiate = (0.11 ± 0.02) + (0.47 ± 0.01) log(r). (4.32) 

At intermediate times L^{t) is growing at a faster rate for longer quenches and 
is not given by a single power-law of the quench rate. 

At late times the texture separation and texture width are well-described by 
the power-law scalings 

Lsep ~ r°-38±°-°\ (4.33) 

~ ^0.47±0.01_ 

The best-fit exponents were derived using equal weighting for all lattice measure- 
ments. As we found no reliable means to estimate the error in individual lattice 
measurements, the errors in the scaling exponents arise from statistical errors 
in the fit. Certainly, power-law scaling is expected from the freeze-out picture. 
However, the Kibble-Zurek mechanism also provides a precise prediction for the 
value of the power-law exponent as derived from equilibrium critical exponents. 
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The model under consideration is a mean field model and the only role of the 
small fluctuations is to seed texture formation. Since the effective mass is con- 
trolled externally and ?^e// ^ tc — t near the critical point, it is expected that 
z/ = i, the mean field value of equilibrium correlation length exponent. This 
expectation was verified by equilibrium lattice simulations. With rj = 1.0 and 
r > 10 the dynamics is overdamped at freeze-out. Therefore /i = 2i/ = 1 and the 
Kibble-Zurek mechanism predicts the scaling 

e/rz(r)~ri, (4.34) 



in obvious disagreement with Eq. ( [4.33|) . 



The larger observed exponent means that for longer quenches topological tex- 
tures are both less numerous and bigger than the Kibble-Zurek mechanism pre- 
dicts. However, slower quenches take longer for the phase transition to complete, 
since as r increases, the order parameter rolls more slowly. In fact, observations 
of the dynamics show that the time to. 9 (measured from the critical point) it 
takes for the order parameter to reach ($^) = 0.9 increases with r in a simple 
way to.9 ~ 1.5r. With more time to evolve between freeze-out and the end of 
the transition, slower quenches allow for more interactions among the texture 
distribution. This suggests that the apparent power-laws observed at the end of 
the transition are not indicative of a single scale, but are a combination of the 
Kibble-Zurek mechanism and the dynamical length scales of the evolving texture 
network. 

Evolution of the power law exponents of the texture length scales is clearly 
indicated in the two lower graphs of Figs. 4.12 and [4.13| . At these early times the 



effect of coarsening dynamics is insignificant and power-law scaling indicative of 
the Kibble-Zurek mechanism is evident, 



Lser, - rO-22±o-oi, (4.35) 
n.i4±o.oi 



Although the power-law exponent of Ly^ is slightly smaller than predicted by the 
Kibble-Zurek mechanism, the difference is not likely to be significant. Textures 
are not well-formed at these early times and the method for measuring Ly^ is 
sensitive to the details of the two-point function. In the middle graphs of each 
figure, the number of textures and the texture width were measured at inter- 
mediate times. Both show the crossover from formation dynamics to coarsening 
dynamics as longer quenches develop a steeper power-law indicative of texture 
interactions. A similar picture of the evolution of Lg^p is given by the average 
number of textures as measured by the total gradient energy. Fig. [4.14| . As be- 



fore, symbols denote lattice measurements and the two solid lines denote best-fit 
power-laws. Early in the phase transition (/) = 0.1 the best fit power-law is 

log{G)eariy = (2.34 ± 0.01) - (0.46 ± 0.01) log(r), (4.36) 
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Figure 4.14: Log-Log plot of the average number of textures determined from the 
gradient energy G = J (Px{d^i ■ d^i) versus quench parameter r for different 

values of the order parameter (p = ($^). 



in agreement with the Kibble-Zurek mechanism. At late times = 0.9, 

hg{G)iate = (3.65 ± 0.10) - (0.94 ± 0.06) log(r). (4.37) 



As in Fig. [4.12| a crossover is apparent at intermediate times = 0.45. The 
late-time power law is slightly steeper than that determined from the topological 
charge because the gradient energy also includes topologically trivial configura- 
tions which decay more rapidly than textures. 



Interactive texture formation 



Figures |4.12| , |4.13| and [4.14| provide clear evidence of the evolution of the power- 
laws characterizing the length scales of the texture distribution. We can use the 
knowledge of the coarsening dynamics of texture interactions to provide a quanti- 
tative understanding of the power-laws observed at the end of the transition. As 



discussed in Sec. ^72| , the phase ordering of a textured system results in a dynami- 



cal length scale characterizing the separation of defects, Li(t) = Qt^. Figure |4?T^ 
shows the dynamics of this length scale for a fast quench (which enters the coars- 
ening regime earlier). The solid line denotes data from lattice measurements. 
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Figure 4.15: Plot of the average number of textures Qit) = J (fx\p{x,t)\ versus 
time for a fast quench r = 1. 



The dashed hne is the best-fit power law of the late-time coarsening dynamics, 

Q{t) ~ ^4 

The coarsening exponent a = —0.72 agrees favorably with previous studies 
which found a = — |. In the instantaneous quenches of ||8^ was identified 
with the initial correlation length. In a dynamical quench it is more natural to 
identify with the freeze-out correlation length. Combining the Kibble-Zurek 
mechanism with phase ordering dynamics we suggest that the evolution of the 
scale characterizing the separation of topological textures is given by 

-'-'Sep it) = reeze reeze )Ht-t freeze)-', (4.39) 

where ^ freeze ~ is the length scale determined by the Kibble-Zurek mechanism, 
tfreeze ~ IS the time from freeze-out and t is the time from the critical point. 
Although Eq. ([4.39|) appears slightly complicated, its form is tightly constrained 



and the only free parameter is a, the exponent of the freeze-out correlation length. 
To test Eq. ([4.39|) we match the expected number of defects Q ~ -p— at time 



t = to.9 ~ 1.5r against the number of defects determined from the numerical 
simulations. A plot of the data and the best-fit length scale is shown in Fig. |4.16. 



Filled triangles denote lattice measurements while the solid line denotes the the- 
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Figure 4.16: Plot of the average number of textures Q = J d^x\p{x) versus quench 
parameter r. 



oretically expected number of textures computed from Eq. ( [4.39|) . To obtain the 
fit one data point is used to normahze the scale and a value of a is determined 
by minimizing x^- By using, in turn, each data point as a normalization a spread 
in the value of a is obtained. The best-fit value of a is taken as the average over 
all normalizations and the error is the standard deviation. As determined from 
these fits, a = 0.24 ± 0.02 in excellent agreement with a = 0.25 expected from 
the Kibble-Zurek mechanism. 

The analysis of the dynamics of L^, the length scale characterizing the size 
of topological textures, is more complicated. Fig. [4.17| provides an example of 
texture width dynamics for a single quench with r = 80. Solid triangles de- 
note lattice measurements and the straight line is the best-fit power law of the 
coarsening dynamics. 



Lw{t — tc) ^ {t — tc 



-0.62 



(4.40) 



The coarsening exponent a = —0.62 is close to the simple scaling exponent 
a = — |. In the simulations the texture width was observed to grow in time with 
a power law near the simple coarsening dynamics predicted for an overdamped 
model L ~ t2. However, this growth rate is different from previously reported 
studies of the late-time dynamics of the texture width. Consistent with previous 
80[, it is possible that Lyj{t) reaches its asymptotic dynamics only at later 



studies 
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Figure 4.17: Log-Log plot of the texture width versus time for quench parameter 
r = 80. Time t is measured from the critical point, tc = 80. 



times. It is also possible that Lyj{t) is contaminated at early times by topologi- 
cally trivial spin configurations whose decay then dominates the dynamics. The 
resolution of this issue requires a more accurate method of counting textures e.g. 
by looking explicitly at the winding around the vacuum manifold. Whatever the 
mechanism, Fig. [4.17| shows that can evolve significantly between the early 
and late periods of the phase transition. We suggest that the evolution is given 

by 

Lw{t) = i freeze + ~ t freeze) ^ • (4-41) 

As before, freeze ~ 'T^ is the length scale determined by the Kibble-Zurek mecha- 
nism, t freeze ~ is the time from "freeze-out" and t is the time from the critical 
point. To test Eq. ( 4.41 ) we match the expected width of the textures at time 



t = 1.5r against the width determined from the numerical simulations. A plot 
of the data and the best-fit length scale is shown in Fig. [4.18| . Filled triangles 
denote lattice measurements. The solid line denotes the theoretically expected 
width computed from Eq. ( [4.41|) . The fits were obtained in the same way as 



for Lsep- As determined from these fits, f3 = 0.22 ± 0.07 also in agreement with 



P = 0.25 expected from the Kibble-Zurek mechanism. It is clear from Figs. |4.1(: 



and [4.18| that the power-laws of Eq. ( |4.33|) can be explained by incorporating 



texture dynamics. Eqns. ( [4.39| ) and (|4.41|) offer a conservative, simple guess to 
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Figure 4.18: Plot of the average texture width versus quench parameter r. 



the proper combination of formation and coarsening dynamics. A more com- 
plete analysis, with better techniques for identifying textures, would fit Lgepif) 
and Lyj{t) to functions with both general formation and coarsening exponents. 
Nonetheless it is remarkable that these simple equations are in good agreement 
with the Kibble-Zurek mechanism. Even if the exact form of Eqns. (|4.39|) and 
( |4.41| ) changes, the observation remains that the length scales characterizing the 
topological defect distribution at the end of the phase transition are a combina- 
tion of defect interaction and formation dynamics. 

It is not surprising that defect interactions can influence the properties of the 
defect network before the end of the transition. Although the long-wavelength 
dynamics is very slow until after the freeze-out time t/ree^e ~ (in the over- 
damped case), the time to complete the phase transition is to. 9 ~ t allowing 
plenty of time for even partially formed defects to influence each other. In pre- 
vious numerical simulations in 1 and 2 dimensions [0, |l^ this effect was not 
noticeable (in overdamped evolutions) because the interactions between defects 
were extremely weak. In a recent simulation of global vortex formation in 3 di- 
mensions 1 19 1 one might expect an observable change in the scaling exponents 
since the defect interactions are much stronger. However, the transition was at 
relatively high temperature and the defects were at least partially screened by 
thermal fluctuations. 
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Texture coarsening and the Kibble-Zurek mechanism 



The close agreement between Eqns. ( [4.39|) and ( [4.41| ) and the length scales ob 



served at the end of the transition provides evidence that the a priori independent 
length scales, Lsep and L^, are formed by the same nonequilibrium dynamics of 
the Kibble-Zurek mechanism. That this should be so is not obvious. In fact, the 
Kibble-Zurek mechanism can be questioned even in models which at late times 
depend only upon a single scale. In it was argued argued that the length 
scales defined by the field correlation length and vortex density, in principle, 
derive from different attributes of the power spectrum. Only under restricted cir- 
cumstances does the defect separation follow the field correlation length and scale 



with the quench rate. The initial analysis of the dissipative quench in Sec. [4.3.1 
may support the contention that strong fluctuations can alter the Kibble-Zurek 
scenario. The external quench, however, occurred in the absence of strong fluctu- 
ations. In this case, as was the case in Chapter 3, the early-time power spectrum 
is dominated by a single peak with momentum k = kmax in the low-momentum 
modes. When the power spectrum is strongly peaked, the freeze-out correlation 
length ^freeze ~ determines a unique scale at early times to which all other 
length scales are connected. 

Late-time dynamics and experiments 

The strong violation of the scaling hypothesis observed in the coarsening of 2 + 1 
dimensional textured systems suggests the possibility of novel, late-time measure- 



ments of early-time critical dynamics. As discussed in Sec. |4.2| , the coarsening of 
a textured system with an instantaneous quench retains "memory" of a length 
scale the correlation length in the disordered phase. In a nonequilibrium phase 
transition it is the freeze-out correlation length, determined by the Kibble-Zurek 
mechanism, and not the correlation length of the initial state, that scars the 
texture distribution at late times. For overdamped coarsening dynamics at late 
times, 

= constant = ^freeze- (4.42) 

The freeze-out correlation length can be measured by a late-time measurement of 
the (average) texture width and texture separation. It is the particular (and not 
well- understood) dynamics of texture-texture interactions that produces strong 
scaling violations and unusual memory effects. In contrast, experimental probes 
of vortex formation in nonequilibrium phase transitions of ^He and ^He are ham- 
pered by the details of the interacting vortex distribution. Because these systems 
approach a scaling solution at late times, vortex interactions work to erase any 
memory of the initial state. This is a particular problem in the ^ife experiments 
where the vortex network is experimentally observable only at later times and 
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the number of initial defects must be inferred using detailed assumptions about 
the decay of the vortex tangle. 

Textured systems in 2 + 1 dimensions can be created in the laboratory. If their 
dynamics are approximately described by Eq. ( 4.19 ) then late-time experimental 
measurements of the Kibble-Zurek mechanism are possible. Topological textures 
can also be formed during phase transitions in 3 + 1 dimensions in superfluid ^He 
and the early universe. In 3 + 1 dimensions without additional higher derivative 
terms in the action, topological textures are unstable. In this case the dynamics of 
coarsening is complicated by texture collapse. However, as long as textures exist 
in the system scaling violations and associated dynamical memory are possible. 



4.5 Summary 

In this chapter, we exploited the multiple length scales of a textured system to 
study a nonequilibrium phase transition of an overdamped classical 0(3) model 
in 2 + 1 dimensions. We introduced a dissipative quench model in which the 
dynamics of the phase transition occur through the dissipative cooling of the 
scalar fields, initially at high temperature. Our analysis of the topological tex- 
tures formed in the dissipative quench seems to indicate a departure from the 
Kibble-Zurek mechanism. Unfortunately an understanding of these results is 
complicated by the presence of large thermal fluctuations and is not yet com- 
plete. We then introduced an external quench model in which the dynamics of 
the phase transition proceed through the controlled change of a time-dependent 
effective mass. 

In the external quench at very early times, we identified power-law scaling 
characteristic of the Kibble-Zurek mechanism and a single freeze-out scale. This 
suggests that the multiple length scales characteristic of the late-time ordering 
of a textured system derive from the critical dynamics of a single nonequilibrium 
correlation length. When observed near the end of the phase transition, we found 
the scaling of the texture separation Lsepij) and the texture width L^(r) result 
instead from a competition between the length scale determined at freeze-out and 
the ordering dynamics of a textured system. We expect that this observation 
is not restricted to systems of topological textures or to low dimensions but is 
relevant anytime defect interactions are significant. It is not surprising that defect 
interactions can significantly modify the defect distribution even before the end 
of the phase transition. Well before defects are fully formed, they frustrate the 
system from its true minimum energy ground state and interact with each other. 

The Kibble-Zurek mechanism provides a useful connection between critical 
dynamics and the length scales of the topological defect distribution observed at 
the end of the phase transition. However, power-law scaling provides a rough and 
rather opaque window into the complicated nonequilibrium processes of the phase 
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transition. As demonstrated explicitly here, power-law scaling can also arise from 
completely different mechanisms such as defect interactions. To understand the 
characteristics of the topological defect distribution at the end of the transition 
it is therefore important to look closely both at the nonequilibrium dynamics of 
the phase transition and the evolution of the defect network. 
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Chapter 5 



Black hole phase transitions 



5.1 Introduction 

As discussed extensively in previous chapters, thermal fluctuations can induce 
phase transitions in which the zero-temperature degrees of freedom are reorga- 
nized into a qualitatively different form. While there are many familiar examples 
such as the boiling of water, phase transitions also occur in more exotic sys- 
tems like spacetime geometry. In this chapter, which represents work in progress 



26| , |27| , pq) , we study a spacetime phase transition evident through the sponta- 
neous formation of a black hole in equilibrium with a thermal environment. We 
seek the answer to two important questions: 

1. Why does a black hole phase transition occur? 

2. By what dynamical scenario does the phase transition take place? 

To answer these questions we will necessarily cast a wide net, drawing insight 
from diverse subjects ranging from Euclidean quantum gravity and string theory 
to vortices in condensed matter. 

The abstract nature of the black hole phase transition is a departure from 
the more physical systems considered in earlier chapters. No laboratory is yet 
equipped with the tools necessary to probe the formation of a black hole in a 
phase transition of thermal spacetime, although such situations may have existed 
in the very early universe. Here, our motivation is not the modeling of an ex- 
isting physical system. Rather we intend to use our knowledge gained from the 
previous study of quantum and classical phase transitions to peer into the compli- 
cated workings of general relativity and semiclassical gravity. General relativity 
is a highly nonlinear theory and comparatively little is known about its detailed 
behavior away from exact solutions. For example, in the extremely energetic en- 
vironment of the early universe, general relativity may exhibit a disordered phase 
dominated by black holes, wormholes, geons and other nonperturbative gravita- 
tional excitations. At even higher energies near the Planck scale it has long been 



speculated that spacetime appears as a foam or froth [Q. While the picture of 
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spacetime foam is vivid, it is very hard to realize in quantitative detail: in part 
because a theory of quantum gravity does not yet exist. However, it is our hope 
that understanding the phase structure of semiclassical gravity provides at least 
a preliminary step on the road to quantum gravity. 

Black holes bear some similarity to topological defects. They are both stable, 
nonperturbative solutions of the classical theory. Both black holes and topologi- 
cal defects carry with them a remnant of the high temperature phase; symmetry 
is restored in the core of topological defects and we are likely to find a quantum 
phase of spacetime in the high-curvature region near the black hole singular- 
ity. In addition, in results reported in this chapter, we show that the black 
hole phase transition is qualitatively similar to the defect-mediated Kosterlitz- 
Thouless phase transition in condensed matter and the Hagedorn transition in 
string systems. 

The analogy between black holes and topological defects is only part of a 
larger relationship existing between condensed matter physics and "fundamen- 



tal" physics such as general relativity, particle physics and cosmology [0, ^ [13 
Such apparently disparate subjects are unified by the complex organizing behav- 
ior of their constituent elements, whether atoms and molecules in condensed 
matter, or spacetime itself in quantum gravity. Condensed matter systems are 
usually more completely understood than their high energy counterparts and we 
can exploit these analogies to illuminate the behavior of systems which are oth- 
erwise experimentally and theoretically intractable. For example, Bose-Einstein 



condensates may provide a testable model of black hole Hawking radiation |Q . 
The experimental tests of physics at high energy scales are (and are likely to 
remain) relatively few. In this environment, the analogies between condensed 
matter and other physical systems are increasingly important in our understand- 
ing of traditional fundamental physics such as semiclassical gravity. 

Organization 

The organization of this chapter is as follows. In Sec. p.2| we review the thermo- 



dynamics of the black hole phase transition, focusing on the calculation of the 
semiclassical free energy of a black hole in thermal equilibrium. In Sec. we 



review an atomic model of the black hole, first introduced as a quantum model 
of black hole microstates. In original work we use this model to provide a statis- 
tical mechanics of the phase transition and highlight the important role of black 
hole entropy. In Sec. |5.4| we discuss the nonequilibrium dynamics of the black 
hole phase transition. We highlight the inadequacy of the homogeneous nucle- 
ation theory of first-order phase transitions and explore examples from condensed 
matter and string theory which appear qualitatively similar to the black hole sys- 
tem. A summary is provided in Sec. |5.5| . In this chapter, unless otherwise noted, 
we use Planck units for which h = G = c = 1. 
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5.2 Black hole free energy 



For normal matter, gravitational interactions are universally attractive and a 
self-gravitating system is fundamentally unstable to collapse. For example, a 
nonrelativistic ideal homogeneous fluid with density p and sound speed Vs is un- 
stable to long wavelength density perturbations. Perturbations with wave vector 




kj<J-f (5.1) 



will grow exponentially. This is the Jeans instability. Since gravity cannot be 
screened, gravitational instabilities remain for a system in thermal equilibrium. 
Compress an ideal, isothermal, self-gravitating gas below a critical volume and 
the gas will collapse. 

The combination of quantum theory and general relativity adds additional 
instabilities to a thermal gravitational system. In particular, hot, fiat space is 
unstable to the (quantum) nucleation of black holes [|^. Using the techniques of 
Euclidean quantum gravity, the nucleation of black holes was identified through 
the discovery of a Schwarzschild instanton contributing an imaginary piece to the 
free energy of the system. The nucleation rate is maximum for a black hole with 
mass M = and is (approximately) given by 



A black hole nucleated with temperature T is in unstable equilibrium with a 
thermal environment of the same temperature. If, in a fluctuation, the black hole 
absorbs a small amount of radiation, its Hawking temperature decreases (black 
holes generally have negative speciflc heat). As the black hole grows it becomes 
colder still, absorbing more radiation and eventually engulflng the system. Thus, 
although r is small except near the Planck scale, the negative speciflc heat of 
nucleated black holes renders the canonical ensemble of hot, flat space ill-deflned. 

Black hole systems become thermodynamically stable with the addition of 
special boundary conditions or in spacetimes with a negative cosmological con- 
stant [0. In the following we flx the temperature T on an isothermal boundary 
of radius r containing a black hole of mass M. In equilibrium, the Hawking 
temperature measured on the boundary must equal the boundary temperature, 

T(r) = \^ , =. (5.3) 



Eq. (|5.3|) admits two real, nonzero solutions for the mass: a smaller, unstable 
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black hole with mass Mi and a larger, stable black hole with mass M2, 



Mo 



SttT 
r 
2 
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The isothermal boundary renders M2 thermodynamically stable because of the 
temperature redshift. A fluctuation that increases M2 also increases the temper- 
ature of the black hole as measured on the boundary, giving M2 a positive specific 
heat. Surprisingly, for 
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(5.6) 



no real value M can solve Eq. ( |5.3| ) and no black hole can exist in the box. It 
therefore appears that as the temperature on the boundary is increased from 
T = 0, a phase transition to a black hole spacetime occurs at T = T^. 

To further elucidate the thermodynamics of the black hole system and the 
nature of any potential phase transitions we consider the canonical partition 
function defined through an Euclidean path integral p6[ . 



Z[P] = / D[g]e 



_IeM. 



(5.7) 



where we have temporarily restored the h dependence in anticipation of the semi- 
classical limit. The Euclidean action is obtained from the Lorentzian Einstein- 
Hilbert action (with boundary terms) through the Wick rotation, t —ir. The 
functional integration is taken over real Euclidean metrics, periodic in imaginary 
time coordinate r with period equal to the inverse temperature (3. Apart from 
artificial toy models, the full functional integral is intractable. However, in the 
semiclassical limit (/i — 0), the integrand is highly peaked around metrics that 
minimize the classical Euclidean action. In the semiclassical limit we evaluate the 
partition function for the system consisting of a single black hole in a finite cavity 
with boundary topology x S"^ (the partition function beyond the semiclassical 
approximation was considered in [0). The action is 



(5.^ 



where K is the trace of the extrinsic curvature of the boundary and 7 is the 
determinant of the induced three-metric. The free energy of a single Schwarzschild 
black hole of mass M within the cavity is given by F = I3~^Ie where the Euclidean 
action Eq. ([5.8[) is evaluated for the metric 



ds^ 



1 



2M 



2M 



(5.9) 
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Figure 5.1: Plot of the free energy vs. mass of the black hole system. 



and r is a Euclidean time coordinate with period (3. Normalized so that F 
when M = 0, the semiclassical free energy is 



F(M, r, T) = r -rJl - — -47r M^T. (5.10) 



In Fig. |5.1| we plot the free energy F{M) at various temperatures for a system 
with box size r = 10. From the top down, the first curve is for temperature 
T = 0.01, below the critical temperature Tc = 0.021. The next curve is for T = T^. 
The lower three curves are for T = 0.25, T = 0.4 and T = 0.5 respectively. 
For temperatures below Tc no black hole is present. At T = Tc an extremum 
appears at M = Mi = M2 (where the free energy is flat). For temperatures 
above there are two extrema, the unstable black hole with mass Mi and the 
stable black hole with mass M2, in agreement with the simple arguments at the 
beginning of the section. Figure ^Tl] also suggests that the transition to a black 
hole spacetime above Tc occurs discontinuously. Above Tc a finite mass black 
hole can be nucleated in equilibrium with the walls of the box. 

The nucleation of a stable black hole at T = does not necessarily signal 
a phase transition. A system in thermodynamic equilibrium always resides in a 
state of lowest free energy. A phase transition occurs only if the free energy of 
the system with the black hole is lower than the free energy without the black 
hole. At temperatures below Tc no black hole is present and the free energy is 
approximately that of a box filled with thermal gravitons (hot, flat space), 

Fhfs ~ -TV. (5.11) 

The free energy of hot flat space is negative. Therefore a phase transition from 
hot flat spacetime to a black hole spacetime can only occur at temperatures 
T > Tc for which F{M2) < F^fg. In Fig. ^]1] the temperature is low enough that 
-^/i/s ~ and a black hole phase transition occurs when F{M2) < 0. 
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The nucleation of a black hole in equilibrium with the walls of the box does not 
appear to be the result of classical gravitational collapse. If the box is uniformly 
filled with massless thermal radiation, we approximate the collapse temperature 
as the temperature for which the Schwarzschild radius for the thermal energy of 
the radiation is equal to the box size, 

TLiapser' r. (5.12) 

Thus Tcoiiapse ~ qualitatively distinct from the nucleation temperature. 

The phase transition does not occur at a high temperature characteristic of 
quantum gravity. In fact, as the size of the box is increased the critical tempera- 
ture decreases and is arbitrarily low for arbitrarily large boxes. The stable black 
hole formed at temperatures above Tc is large, with a mass M2 on the order of 
the size of the box. These considerations seem to be paradoxical. How can such 
a large energy [E ~ M) fluctuation actually lower the free energy of the system 
at such low temperatures T <C M? The answer lies in the enormous entropy 
black holes hold within their horizon. The free energy results from a competition 
between the internal energy and the entropy, F = E — TS. For a black hole, 
E = M and the entropy is proportional to the area of the horizon, S = AttM'^. 
Because the entropy is growing with mass faster than the energy, there is always 
a critical temperature above which the entropy completely compensates for the 
energy cost of making a black hole and the black hole spacetime is the lowest free 
energy state. 



5.3 Black hole atoms 

The black hole phase transition is entropically driven. Therefore, an understand- 
ing of the nature of black hole entropy offers potential insight into the details of 
the transition. In this section we study a toy model for black hole microstates, 
focusing on their implications for the thermal black hole system. 

The origin of black hole entropy has been an outstanding problem since Beken- 



stein first introduced the concept |98|. A complete resolution likely requires a 
consistent theory of quantum gravity, which has so far proved elusive. However, 
just as semiclassical reasoning such as the Bohr model was important in the early 
development of quantum theory and the interpretation of atomic spectra, a sim- 
ilar approach may be fruitful in understanding some of the microscopic features 



of black hole entropy |99]. It is not our intention to survey the large number of 



semiclassical black hole models. However common to many is the quantization 



of the horizon area into equally spaced levels (see ||100|| and references therein), 

A = an; n = 1,2 . . . (5.13) 



81 



where a is dimensionless but as yet unspecified and the area is given in units of 
Planck area. For a black hole of mass M the entropy, 



5=4= 47rM2, 
4 

and the mass of the black hole is also quantized, 

M = 7Vn, 



(5.14) 



(5.15) 



where 7 = -^Z^- If the energy levels have degeneracy g{n), the black hole 
entropy may simply count the number of available microstates for a fixed energy 
level n, S = In g{n). In this case. 



9[n) 



Since g{n) must also be an integer, a is restricted. 



a 



4\nk: 



k = 2A... 



(5.16) 



(5.17) 



Equations ( |5.15| ), ( |5.16|) and ( |5.17| ) define a semiclassical black hole model in 
analogy with the Bohr model of an atomic system. If = 2 the large degeneracy 
g{n) can be thought of as the number of ways to make a black hole in the nth 
level by starting in the ground state ||101||, though there are other interpretations 



|102| , |103|| . Our interest in the atomic black hole model is its behavior in thermal 
equilibrium. 

The partition function for the quantum black hole atom in the canonical 
ensemble is 



T 



(5.18) 



n=0 



Without modification, the sum in Eq. ( p.l8| ) does not converge. Convergence is 
obtained upon analytic continuation but the partition function acquires a imagi- 
nary piece | 103| |. In light of the discussion of the previous section this is not at all 
surprising. Without either special boundary conditions or special spacetimes, the 
thermodynamics of black holes is not well-defined. In fact, Im(Z) is due precisely 
to the nucleation of black holes. To obtain a well-defined and real partition func- 
tion we take a new approach and place the quantum black hole atom into a box 
of radius r. The box is realized as a sharp cutoff in the energy levels accessible 
to the black hole. 



472 



(5.19) 
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Figure 5.2: Plot of the average energy E vs. temperature T for the black hole 
atom. 
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Figure 5.3: Plot of the specific heat vs. temperature T for the black hole 
atom. 



This is reasonable as the box functions to remove energy levels with Schwarzschild 
radius larger than r. The canonical partition function for the quantum black hole 
atom in a box is 



47^ 



7^" 
T 



(5.20) 



n=0 



In Figs. [5.2| , [5.3| and |5.4| we plot the average energy, specific heat and entropy as 
a function of temperature for a quantum black hole atom in a box with radius 
r = 10. The plots were made with rescaled variables T — >■ 7T, r — > 7r and 
for the particular choice k = 2. All thermodynamic quantities were calculated 
using the partition function Eq. (|5.20|) . The inspection of Figs. ^.2[ - ^.4| reveals 
a sharp transition from the ground state to a highly excited state in the black 
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Figure 5.4: Plot of the entropy S vs. temperature T for the black hole atom. 



hole atomic system. This is reminiscent of the phase transition behavior that 
we saw in the previous section. To quantify this behavior consider the effective 
Boltzmann factor of level n, 



(5.21) 



/(n) = — + nlnk. 



For small n, f{n) is negative and higher energy levels are suppressed, as is usually 
the case. However, since the degeneracy grows as ra, there is always an energy level 
nc(T) defined by /(ric) = beyond which degeneracy compensates for the higher 
energy. Levels beyond are enhanced, not suppressed. At low T, u^iT) > n^ax 
and these enhanced levels are not part of the allowed spectrum. As T increases 
nc{T) decreases and a transition to the higher levels occurs when nc{T) = Umax- 
If we adopt this picture of the black hole phase transition the critical temperature 
is 



Tr 



1 



27rr 



The average energy at the critical point is 



~ e 



/K{Tc)) 



(5.22) 



(5.23) 



These equations reveal two important details. First, although the numerical 
factors are slightly different they have the same qualitative structure as their 
thermodynamic counterparts, Eqns. ( |5.6| ) and ( |5.4| ). Second, neither equation 
contains the one free parameter k of the black hole atom. In fact Eqns. ( |5.22|) 
and ( 5.23| ) are largely independent of the details of the spectrum of the black hole 
atom and apply even when the horizon area is quantized into nonuniform levels. 

At first glance, our study of the quantum black hole atom in thermal equi- 
librium appears to only slightly advance our understanding of the black hole 
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phase transition. We have simply reaffirmed the fact that a system with finite 
energy (hmited by the boundary conditions) and obeying a thermodynamic rela- 
tion S = AttE'^ will always have a transition to the highest allowed energy state, 
the details of which are largely model independent. But in principle the quantum 
black hole atom allows us to go farther. We can now imagine treating dynamical 
processes such as emission and absorption much as we do with more common 
atomic systems. In fact only two levels, the ground state and the highest allowed 
excited state, are significantly populated below and above the transition and we 
may further approximate our system as a quantum two-level atom, a popular sys- 
tem of study. Since the quantum black hole transition occurs as the excitation 
of a high energy state, the dynamics of this atomic excitation approximate the 
dynamics of the black hole phase transition. We hope to report on further work 
on this topic in the near future p6 |. 

Our analysis of the quantum black hole atom in thermal equilibrium under- 
scores the entropic nature of the black hole phase transition. This result is useful 
in itself. Following the dynamics of the black hole phase transition through a 
fully nonequilibrium formulation of semiclassical gravity is a very hard problem 
to which, at the moment, there is no direct method of attack. However, the study 
of simpler systems with a similar entropic transition is likely to yield insight into 
the dynamics of the black hole phase transition. In the next section we turn our 
attention to these dynamical issues. 



5.4 Nonequilibrium dynamics 



Figure appears qualitatively similar to the free energy of a system undergoing 
a (strongly) first-order phase transition. The field theoretic treatment of the dy- 
namics of first-order phase transitions is based upon the homogeneous nucleation 
theory developed by Langer (see for example ||104|| ). In homogeneous nucleation. 



widely separated spherical bubbles of the stable phase nucleate in a background 
of the unstable phase. If the bubbles are larger than a critical radius, the volume 
energy of the stable phase inside the bubble is less than the surface energy and 
the droplet will grow. The phase transition completes as droplets of the stable 
phase expand to fill the volume of the system. In this picture, the black hole 
phase transition occurs when an unstable black hole with mass Mi nucleates 
(with probability P ~ g-^^^^-^i)) and grows to form the stable mass M2 through 
the absorption of thermal radiation. There are reasons to believe, however, that 
homogeneous nucleation is not the correct description, as we indicate below. 

The free energy of the black hole system is calculated under the assumption 
of spherical symmetry, adequate only for a single black hole. In equilibrium, a 
single black hole is preferred because it maximizes the entropy. For example, in 
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a state with two black holes, 



Sma + 'S'm, r^ml + ml < Sraa+m^ ~ ("^a + "^f,)^- (5.24) 

However the dynamics of the phase transition may involve multiple black holes. 
Black holes are strongly interacting gravitational systems, unscreened by thermal 
fluctuations. It is possible that the dynamics of the phase transition proceeds 
by the (exponentially more probable) nucleation of small black holes with mass 
m -C Ml which then merge to form larger holes. Phase transitions in which there 
are strong interactions between bubbles of the new phase lie beyond the scope of 
homogeneous nucleation. 



5.4.1 Black hole gas 

To study the possibility that the dynamics of the phase transition proceeds 
through the nucleation and merger of small black holes we propose a model in 
which black holes are treated as a gas of gravitationally interacting particles. The 
number of particles is not flxed, but changes through the stochastic nucleation of 
small black holes, mergers and Hawking evaporation. To incorporate the inter- 
action of black holes with the thermal environment the mass of each particle is 
not constant but changes in proportion to the free energy of a single black hole. 



dm{t) SF 
dt 6m 

From F{m) given by Eq. ( ^.10|) we obtain the phenomenological relation 

dm{t) 1 



(5.25) 



dt 



2m 



+ 87rmT, (5.26) 



where T is the temperature and r is the size of an effective boundary. With 
a time-dependent mass given by Eq. ( |5.26| ) black holes are stable only above a 
critical temperature, as we expect from previous discussions. The difference and 
advantage of this model is that above the critical temperature, a stable black 
hole can form by small mergers, in addition to a single large fluctuation. A 
conceptually similar approach but with very different emphasis was considered 
in ||105|| . We hope to report on the analysis of this black hole gas model in the 



near future . 



5.4.2 Entropic transitions 

Another approach to the dynamics of the black hole phase transition is through 
the study of similar physical systems. The black hole phase transition is driven 
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by the large amount of black hole entropy S ^ E"^. A thermodynamic relation 
where the entropy grows rapidly with energy is very unusual. For example, 
a classical ideal gas has entropy 5* ~ ln{E). However there are examples of 
condensed matter systems that exhibit similar entropic transitions. In particular 
we examine the Kosterlitz-Thouless (KT) transition in a global 0(2) model in 



2 + 1 dimensions B and the Hagedorn transition [106] in string systems. 



Kosterlitz-Thouless transition 

We consider a global 0(2) scalar field model with Hamiltonian, 

H = ld'x Q|V0? + ^(0^ - ■ (5.27) 

This model admits vortex topological defects. The energy of a vortex of single 
winding is 

E ^ nr]^ (^n ^ + Ar/^a^^ , (5.28) 

where a ~ is the vortex size and R is the size of the system. If i? ^ a the 
energy of the vortex is dominated by gradient energy, the first term in Eq. ( p. 281 ) . 
If a vortex can be nucleated anywhere in the system then the entropy 

S = \n(^^^ , (5.29) 

and the free energy of the system with a single vortex is 

F = {nr]^ -2T)\n(^^y (5.30) 

Inspection of the free energy reveals that above the critical temperature 

Tc=^, (5.31) 

it is thermodynamically favorable to nucleate vortices. This is the KT transition, 
a transition from (algebraic) order to disorder. The dynamics of the KT transi- 
tion proceeds through the interaction of vortices. At low temperatures vortices 
exist as a low density of tightly bound vortex-antivortex pairs with zero net topo- 
logical charge. As the temperature increases, both the density of pairs and the 
average vortex-antivortex separation also increase. In addition, as the separation 
increases, thermal fluctuations are more effective in screening the interaction 
between defects in a vortex-antivortex pair. When the critical temperature is 
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reached, screening is complete and the vortex-antivortex pairs unbind leaving a 
gas of essentially free topological defects. 

The KT transition exhibits both similarities and differences when compared 
with the black hole system. Analogous to the black hole transition, the nucleation 
of a vortex with energy E produces a large amount of entropy S ^ E which drives 
the phase transition. In distinction however, the vortex entropy results from the 
possible locations of the vortex and not, as in the black hole, from internal states. 
The natural analogy of an antivortex in the black hole system is a white hole, a 
spacetime region that expels all worldlines. A white hole is the time-reverse of a 
black hole, just as (e.g. in superfluid helium) an antivortex is the time-reverse of 
a vortex. The analogy of the KT transition suggests that the black hole phase 
transition might be understood as the unbinding of black hole-white hole pairs. 
It is far from clear that this viewpoint is correct. However, at the very least, the 
analogy of the KT transition suggests a new angle towards the dynamics of the 
black hole phase transition. 



Hagedorn transition 

A system of fundamental closed bosonic strings in thermodynamic equilibrium 



possesses an exponential density of states ||106 | 



z/(E) ~exp(/?H^), (5.32) 

where jSn is the model and dimension dependent Hagedorn temperature. Above 
the Hagedorn temperature the exponential density of states renders the partition 
function ill-defined and a physical description of the system is unclear. 

Insight into the interpretation of the Hagedorn transition is obtained by ex- 
amining the relationship between string theory and Quantum Chromodynamics 
(QCD). QCD is the (nonabelian) SU(3) gauge theory of strong interactions. The 
fundamental degrees of freedom are three colored fermionic quarks interacting 
via bosonic gluons. QCD displays asymptotic freedom, at high energies the QCD 
coupling is weak and quarks are effectively free. However, at low energies QCD 
is strongly coupled and quarks exist only in tightly bound color singlet states, 
hadrons and mesons. Between the two regimes, QCD is expected to undergo 
a deconfinement phase transition at which hadrons and mesons melt into their 



constituent quarks. The large-N ||107|| and strong coupling ||108|| limits of QCD 
may be described by a fundamental string theory and the link between QCD and 
string theory provides a physical picture of the Hagedorn transition: the prolifer- 
ation of string states at the Hagedorn temperature corresponds to the emergence 



of quark color degrees of freedom at deconfinement ||109| , 110 , 111 . 

Surprisingly, string theory also offers the possibility that the Hagedorn tran- 
sition in large-N SU(N) gauge theory is linked through a hypothesized AdS/CFT 
duality to the formation of a black hole in anti-deSitter (AdS) spacetime (see 
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112|| for a comprehensive review). The thermodynamics of a semiclassical black 



hole system in a spacetime with negative cosmological constant is qualitatively 
similar to the black hole in a box discussed previously in Sec. |^. In particular, a 
black hole phase transition, the Hawking- Page transition [95] occurs at a critical 
temperature 

Te ~ (5-33) 

where A is the (negative) cosmological constant. Exploiting the AdS/CFT du- 
ality, Witten argued |[L13|| that the enormous entropy released in the deconfine- 
ment transition of a supersymmetric gauge theory defined on the boundary of 
AdS spacetime has a dual bulk interpretation as a black hole forming through 
the Hawking-Page phase transition. 

A Hagedorn transition also occurs in systems of cosmic strings and condensed 
matter vortices. In this case, the exponentially growing density of states simply 
reflects the "wiggles" present in each string and the Hagedorn transition signals 
the abrupt formation of infinite string. The thermodynamic behavior of strings 
in a classical U(l) scalar field theory with symmetry breaking was examined in 
10, 0. The formation of infinite string was observed at the field critical point, 
suggesting that, as in the KT transition, strings play a fundamental role in the 
phase transition. 

Driven by the exponential density of states, the Hagedorn transition in a 
classical U(l) field theory provides a tractable analogy to the black hole phase 
transition. Visible as the formation of infinite string, we propose to study the 
detailed dynamics of this transition |2^. A U(l) scalar theory provides a simple 



realization of a system exhibiting vortex excitations and also belongs to the same 
important universality class as '^He and type-II superconductors. The use of a 
classical field theory allows for detailed numerical observation of the formation 
and interaction of the nonperturbative string structures important for the transi- 
tion. As discussed above, a Hagedorn transition also describes the deconfinement 
phase transition in large-N or strongly-coupled QCD. Through the AdS/CFT 
correspondence the deconfinement phase transition corresponds to the Hawking- 
Page black hole phase transition. It is therefore possible that the dynamics of 
the formation of infinite string is similar to the dynamics of black hole formation 
in the Hawking-Page transition. Of course, it is impossible to proceed rigorously 
in such a long chain of analogies. However, the dynamics of the formation of 
infinite string in a classical system is an interesting physical problem in its own 
right which, to our knowledge, has not yet been addressed. 



5.5 Summary 

In this chapter we have built the foundation for a thorough analysis of black 
hole phase transitions in semiclassical gravity. We reviewed the semiclassical 
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thermodynamics of the black hole system and determined that the black hole 
phase transition is entropically driven: it occurs because the large entropy of a 
black hole compensates for the energy cost of formation and lowers the total free 
energy of the system. This view was reinforced through the analysis of an atomic 
model of a quantum black hole in thermal equilibrium. In this model, the phase 
transition occurs as the abrupt excitation of a high energy state above the critical 
temperature. Thus, the study of atomic emmision and absorption may provide 
intuition about the black hole system, a direction we are currently pursuing. The 
detailed dynamics of the black hole phase transition remain an open problem. 
We argued that homogeneous nucleation does not apply to phase transitions in 
which there are strong interactions among bubbles of new phase, as is the case 
for the black hole system. We then appealed to the study of similar entropic 
transitions and observed that the black hole phase transition contains elements 
of both the Kosterlitz-Thouless and Hagedorn phase transitions. Reasoning by 
analogy, we argued that the dynamics of the black hole phase transition might be 
similiar to the dynamical formation of long string in the nonequilibrium quench 
of a classical U(l) scalar field theory. Work in these directions is in progress. 
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Chapter 6 



Conclusion 



In this dissertation we have analyzed the nonequihbrium formation of topological 
defects during a second-order phase transition of both a quantum field in the 
early universe and a classical field appropriate for the description of a condensed 
matter system. We also discussed the existence and dynamics of a black hole 
phase transition in semiclassical gravity. 

We first studied the symmetry-breaking phase transition of a quantum scalar 
field in the early universe. Using equations of motion derived from the two- 
loop 2PI-CTP effective action we identified domains in the infrared portion of 
the momentum-space power spectrum. We found that the domain size scales as 
a power-law with the expansion rate of the universe. The observed power-law 
scaling, for both overdamped and underdamped evolution, is in good agreement 
with the predictions of the Kibble-Zurek mechanism. 

We then studied the formation and interaction of topological textures in a 
nonequihbrium phase transition of a classical 0(3) scalar field theory in 2 + 1 
dimensions. We presented two models of the nonequihbrium phase transition, 
a dissipative quench and an external quench. In the external quench and at 
early times, we observed power-law scaling of the length scales of the texture 
distribution in good agreement with the Kibble-Zurek mechanism. However, 
by the end of the transition the length scales of the texture distribution result 
from a competition between the length scale determined at freeze-out and the 
ordering dynamics of a textured system. Therefore a quantitative understanding 
of the defect network at the end of the phase transition generally requires an 
understanding of both critical dynamics and the interactions among topological 
defects. 

Finally, we studied a black hole phase transition in semiclassical gravity. We 
reviewed the thermodynamics of the black hole phase transition, focusing on the 
calculation of the semiclassical free energy of a black hole in thermal equilibrium. 
We determined that the black hole phase transition is entropically driven. We 
applied a quantum atomic model to the black hole equilibrium system and showed 
that the phase transition occurs as the excitation of a high energy state. We 
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discussed the nonequilibrium dynamics of the black hole phase transition and 
explored similar examples from the Kosterlitz-Thouless transition in condensed 
matter to the Hagedorn transition in string theory. 

The power-law scaling of the defect density, observed here in both quantum 
and classical nonequilibrium phase transitions, demonstrates, at least at early 
times, the viability of the Kibble-Zurek mechanism of topological defect forma- 
tion. However, as mentioned in Chapters 3 and 4, the power-law scaling of the 
defect density is only an indirect and partial verification of the ideas presented 
in the freeze-out scenario. In fact, the exact details of the dynamical process 
through which the field correlation length controls the defect separation remain 
somewhat obscure. Partly this is due to the complicated relationship between the 
field and defect degrees of freedom. To remedy this, one possibility is to seek a 
dual topological defect description of the entire nonequilibrium phase transition. 
At least for certain systems we may be able to phrase the freeze-out scenario solely 
in terms of equilibrium defect distributions and equilibrium relaxation times for 
these defect distributions. Indeed, some tentative steps have already been made 
in this direction ||114| . 



In the context of nonequilibrium phase transitions in quantum field theory, an 
important issue is the development of techniques that would allow for an analysis 
of the entire phase transition. As discussed in Chapter 3, although the two- loop 
approximation of the 2PI-CTP effective action is adequate at early times in the 
phase transition, when long wavelength modes are sampling the spinodal instabil- 
ity, the approximation is too simplistic to connect to late-time oscillations around 
the true vacua and the formation and evolution of topological defects. To ana- 
lyze the entirety of the phase transition in nonequilibrium quantum field theory 
we propose a two-stage model of the phase transition. In the first stage quan- 
tum field fluctuations can be evolved at early times using the known techniques 
of nonequilibrium quantum field theory. In the second stage, when standard 
techniques fail, we can attempt to match the field fluctuations to an effective 
(stochastic) classical Landau- Ginzburg equation. The late-time evolution of the 
Landau-Ginzburg equation allows the field to settle to the true vacuum. A sim- 
ilar program was initiated in Ref. but with neither the focus on topological 
defects nor an explicit solution of the late-time dynamics. The advantage of this 
approach is a picture of the transition that extends from the phase dominated by 
quantum fluctuations to the phase dominated by topological defects. In quantum 
field theory, such a complete picture has never been given. Along these lines we 
continue our research. 
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